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ABSTRACT 

The relativistic wind of obliquely-rotating pulsars consists of toroidal stripes of opposite magnetic 
field polarity, separated by current sheets of hot plasma. By means of two- and three-dimensional 
particle-in-cell simulations, we investigate particle acceleration and magnetic field dissipation at the 
termination shock of a relativistic striped wind. At the shock, the flow compresses and the alternating 
fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength A or the wind 
magnetization a (in the regime cr ^ f of magnetically-dominated flows), shock-driven reconnection 
transfers all the magnetic energy of alternating fields to the particles, whose average Lorentz factor 
increases by a factor of a with respect to the pre-shock value. The shape of the post-shock spectrum 
depends primarily on the ratio X/{ri^a), where r^ is the relativistic Larmor radius in the wind. The 
spectrum becomes broader as the value of X/{r]^a) increases, passing from a relativistic Maxwellian to 
a flat power-law tail with slope around —1.5, populated by particles accelerated by the reconnection 
electric field. Close to the equatorial plane of the wind, where the stripes are symmetric, the highest 
energy particles resulting from magnetic reconnection can escape ahead of the shock, and be injected 
into a Fermi-like acceleration process. In the post-shock spectrum, they populate a power-law tail with 
slope around —2.5, that extends beyond the flat component produced by reconnection. Our study 
suggests that the spectral break between the radio and the optical band in Pulsar Wind Nebulae can 
be a natural consequence of particle acceleration at the termination shock of striped pulsar winds. 
Subject headings: acceleration of particles — galaxies: jets — gamma-ray burst: general — pulsars: 
general — radiation mechanisms: nonthermal — shock waves 



1. INTRODUCTION 

The broadband spectrum of Pulsar Wind Nebulae 
(PWNe), extending from the radio up to the gamma- 
ray band, is usually modeled as synchrotron and in- 
verse Compton radiation from a nonthermal population 
of electron-positron pairs. It is assumed that the emit- 
ting particles are accelerated to a power-law distribution 
at the so-called termination shock, where the momentum 
flux of the relativistic pulsar wind is balanced by the con- 
fining pressure of the nebula. However, it is still unclear 
how the wind termination shock can accelerate electrons 
and positrons (hereafter, simply "electrons" ) to the rela- 
tivistic energies required to power the nebular emission. 

A flat radio spectrum is the common observational fea- 
ture of PWNe, whose spectral flux in the radio band 
scales as F^, ex i/""'', with a^ = 0.0 — 0.3. For the Crab 
Nebula, the best studied example of a PWN, the syn- 
chrotron cooling time of radio-emitting electrons largely 
exceeds the age of the nebula, so one cannot exclude that 
such particles may be "primordial," i.e., acce lerated at 
a very early stage of the PWN evolution ( |Atoyan^^I999j) 
However, observatio ns of wisps in the radio band ( Br ^ 
etenholz et al.|200T ) seem to suggest that radio-emittmg 
electrons are being accelerated now, in the same region 
as those responsible for the optical and X-ray emission. 
In this case, the observed radio spec tral index (a^ ~ 0.25 
for the Crab;'Bietcnholz et al. 1997*) would imply a distri- 
bution of emitting electrons ot the form dN/dE ex E^'p , 
with slope p = 2ar -f 1 ^ 1.5. To explain the radio 
through optical emission of the Crab, the power law of 



shock-accelerated electrons sho uld span at least three 
decades in particle energy (e.g., Lyubarsky 2003). Flat 
electron spectra with slope p < 2 below OeV energies are 
also required t o model the radio em ission of hotspots in 
radio galaxies (Stawarz et al. 20071 and X-ray obse rva- 
tions of luminoias blazar sources ( [Sikora et al. 2009 1 . 

For a power- law particle spectrum with 1 < p < 2, 
most of the particles lie at the low-energy end of the dis- 
tribution, but the plasma energy content is dominated by 
particles at the high-energy cutoff. This contrasts with 
the energy spectrum expected from Fermi acceleration 
in rel ativistic shocks, which normally yields slopes jp > 2 
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(e.g., Achterberg et al. 2001 Keshet fc Waxman|2005 1 , in 
which case low-energy particles dominate both by num- 
ber and by energy. To produce broad spectra with p < 2, 
a mechanism is required that raises the mean electron 
energy, but leaves most of the particles a t relatively low 
energi es. I n the context of pulsar winds, |Hoshino et al.| 
(19921 and Amato & Arons (2006) have proposed that, 
in a wind loaded with ions, resonant absorption of ion 
cyclotron waves by electrons and positrons could gener- 
ate flat distributions downstream from the termination 
shock. However, the required ion injection rate largely 
exceeds what is expec ted from stan dard models of pulsar 
magnetospheres (e.g., Arons|[2007 ) . 

An alternative possibility was discussed by |Lyubarsky| 
( 2003 1 , under the assumption that the flow upstream of 
the termination shock consists of alternating stripes of 
opposite magnetic polarity, separated by current sheets 
of hot plasma (from now on, a "striped wind"). For 
obliquely-rotating pulsars, this is the configuration ex- 
pected around the equatorial plane of the wind, where 
the sign of the toroidal field alternates with the pulsar 



period. If the wind remains dominated by Poynting flux 
till the t ermination shock (which is still an open ques- 
tion, see [Lyubarsky fc Kirk 2001 Kirk &: Skjaeraasen 



20031, shock-compression of the stripes may drive an 



mhilation of the alternating fields, a process known as 
"driven magnetic reconnection." The resulting transfer 
of energy from the fields to the particles could potentially 
generate flat electron distributions, with slope 1 < p < 2. 

The physics of magnetic reconnection can be captured 
self-consistently only by means of multi-dimensional 
particle-in-cell (PIC) simulations. Fully-kinetic PIC sim- 
ulations provide a powerful tool to explore the micro- 
physics of collisionless plasmas, since they can capture 
from first principles the fundamental interplay between 
charged particles and electromagnetic fields. In the con- 
text of relativistic magnetic reconnection in pair plasmas, 
most studies have explored the so-called "undrivcn rccon- 
nection" process (Zenitani fc IIoshino_2001j^2005a^^2007j 
2008 Jaroschek et al. 2UU4), where field annihilation is 
initiated by a transient seed perturbation to an otherwise 
stable current sheet. As discussed above, this is not the 
setup expected at the wind termination shock. Here, it 
is the shock-compression of the flow that steadily drives 
regions of opposite magnetic field polarity toward each 
other, causing reconnection. 

Shock-driven reconnec tion in magnetically-domi nated 
flows has been studied by Petri fc Lyubarsky ( 2007 1 with 
analytical methods and one-dimensional (lU) PIC simu- 
lationsF] They find a condition on the wind magnetiza- 
tion ana the stripe wavelength needed for full dissipation 
of the alternating fields at the termination shock (roughly 
speaking, the post-shock Larmor radius should be larger 
than the separation between current sheets). However, 
in ID all particles gain energy at the same rate, which 
results in a Maxwellian distribution downstream from 
the shock. The generation of nonthermal tails is possi- 
ble only if different particles experience different energy 
gains, which requires multi-dimensional simulations. As 
a preliminary step toward multi-dimensional studies of 
the terminatio n shock in a striped wind, Lyubarsky fc| 
Liverts (20081 have performed 2D simulations of driven 



magnetic reconnection. Two stripes of opposite magnetic 
polarity were compressed by means of an external force, 
which would imitate the effect of a shock. They found 
that driven magnetic reconnection can produce flat non- 
thermal tails, with p ^ I. Yet, their study was limited to 
a single current sheet, and most importantly it did not 
self-consistently account for the presence of a shock. 

In this work, we explore via multi-dimensional PIC 
simulations the acceleration of particles at the termina- 
tion shock of a striped relativistic electron-positron wind. 
The physics of a shock interacting with multiple cur- 
rent sheets is captured self-consistently in 2D and 3D, 
and the transfer of energy from the fleld to the particles 
via shock-driven reconnection is investigated from first 
principles. We find that the alternating fields are com- 
pletely dissipated upon compression by the shock, and 
their energy is transferred to the particles, regardless of 
the prope rties of the fiow. This con trasts with the con- 
clusions of Petri fc Lyubarsky (2007), who found efficient 
field dissipation only tor high magnetizations and short 

^ A similar ID study, but for winds dominated by kinetic (rather 
than Poynting) flux, was performed by Nagata et al. (2008). 



stripe wavelengths. We show that the reduced dimen- 
sionality of their ID model could not correctly capture 
the growth of the so-called "tearing-mode instability," 
which is of paramount importance for field annihilation 
in 2D and 3D configurations. We find that broad parti- 
cle spectra with flat slopes (1 < p < 2) are a common 
by-product of shock-driven reconnection, but the extent 
of the power-law tail depends on the wind magnetization 
and the stripe wavelength. 

This work is organized as follows. In f|2J we discuss the 
setup of our simulations and the magnetic field geometry. 
In fpj the properties of the termination shock are inves- 
tigated for one representative choice of wind parameters. 
We discuss both the shock structure, with focus on the 
physics of shock-driven magnetic reconnection, and the 
process of particle acceleration to nonthermal energies. 
In fpj we study how the particle spectrum downstream 
from the shock depends on the physical conditions in 
the wind, like magnetization and stripe wavelength. Al- 
though we mostly employ 2D simulations, in f|5]we com- 
pare 2D and 3D results, finding good agreement. We 
summarize our findings in S|6] and comment on the appli- 
cation of our results to astrophysical scenarios. 

2. SIMULATION SETUP 

W e use the 3D electromagnetic PIC code TRISTAN- 
MP ( |Buneman||1993[ iSpitkovskylpOOS l to study the ter- 
mination shock of a relativistic striped wind. The shock 
is set up by reflecting a magnetized electron-positron fiow 
from a conducting wall located at a; = (Fig. [I]). The in- 
teraction between the incoming beam (that propagates 
along —x) and the reflected beam triggers the forma- 
tion of a shock, which moves away from the wall along 
+x (Fig. [I]) . The simulation is performed in the "wall" 
frame, where the downstream plasma is at rest. 

To follow the shock evolution for longer times with 
fixed computational resources, we mainly utilize 2D com- 
putational domains in the xy plane, with periodic bound- 
ary conditions in the y direction. In f|5] we compare 2D 
and 3D simulations and show that 2D runs can capture 
most of the relevant physics. For both 2D and 3D do- 
mains, all three components of particle velocities and 
electromagnetic fields are tracked. 

The incoming electron-positron stream is injected 
along —X with bulk Lorentz factor 70. The flow carries 
a strong magnetic field, oriented along the y direction. 
It follows that the shock will be "perpendicular," i.e., 
with magnetic field orthogonal to the shock normal. The 
choice to initialize the field in the simulation plane (i.e., 
oriented along y, as opposed to z) is motivated by the 
agreement between 2D simulations with in-plane fields 
and 3D experiments, as we show in Sj5] The spatial pro- 
file of the injected magnetic field, sketched in Fig. [l| is 



By{x, t) = Bq tanh < - 
= Bq tanh^ 



a + cos 
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and its associated motional electric field is E^ — PoBy. 
Here, /3o is the three- velocity of the injected plasma. This 
particular choice for the magnetic field is motivated by 
analytical studies o f the structure of pulsar winds (e.g., 
Petri fc Kirk 20051). At each given time, eq. (fTl) shows 
that the field profile is periodic with wavelength A (see 




Fig. 1. — Upper panel: poloidal structure of the s triped pulsar 
wind, according to the solution by |Bogovalov| ( |1999| l . The arrows 
denote the pulsar rotational axis (along il, vertical) and magnetic 
axis (along /x, inclined). Within the equatorial wedge bounded by 
the dashed lines, the wind consists of toroidal stripes of alternating 
polarity (see the reversals of Bip), separated by current sheets (dot- 
ted lines). At latitudes higher than the inclination angle between ii 
and /i (i.e., beyond the dashed lines), the field does not alternate. 
Lower panel: simulation geometry. For 2D runs, the simulation 
domain is in the xy plane, oriented as shown in the upper panel 
(so, X = —f and y = —if>). The incoming flow propagates along 
—X, and the shock moves away from the reflecting wall (located 
at X = 0) toward +x. The magnetic field lies in the simulation 
plane along the y direction, and its polarity alternates with wave- 
length A (red stripe for By = +Bq, blue for By = —Bq). A net 
stripe-averaged field {By)x > is set up by choosing red stripes 
wider than the blue stripes. For the pulsar wind sketched above, 
{By)\ > is realized below the equator. 



Fig. [T]). Current sheets, where the magnetic field van- 
ishes, are identified by the condition C ^ 0, and their 
half-thickness is A = X5/2tt ^ A. Across each current 
sheet, the field reverses its sign from +Bo to —Bq (or vice 
versa) , where Bq is the field strength in the region outside 
the current sheets, which we shall call "cold wind." In 
the following, we parameterize the field strength via the 
so-called magnetization parameter a = Bq / ^-K^fomncoc^ , 
taken to be much larger than unity. Here, m is the elec- 
tron (or positron) mass and n^o is the density of parti- 
cles in the cold wind. By defining the relativistic skin 
depth in the cold wind c/wp = y'7omc^/47re^n^, and 
the relativistic Larmor radius r^ = jQinc^ /cBq, we can 
rewrite the magnetization parameter as cr = (c/ajp)^/r|^. 
Finally, the parameter a (—1 < a < 1) is a measure 
of the net field (i.e., averaged over one wavelength A), 
such that {By)x/Bo = oi/[2 — \a\). Although the mag- 
netic field intensity in the cold wind is always i3o, the 
wavelength-averaged field is not necessarily zero, because 
two neighboring stripes generally have different widths 
(see Fig. IT]). In the pulsar wind, one expects {By)x — 
(which corresponds to a = 0) only in the equatorial plane 
(see upper panel in Fig. [Tl. We remark that a net field 
{By)\ 7^ does not play the role of a guide field, which 
in our geometry would correspond to a uniform compo- 



nent of magnetic field in the z direction. No significant 
guide field is expected at the termination shock of pulsar 
winds, and the physics of shock-driven reconnection in 
the presence of a guide field will be discussed elsewhere. 

In the cold wind, each computational cell is initial- 
ized with two electrons and two positrons, with a small 
thermal spread 8c = kTc/mc^ = 10^''. We have also 
performed limited experiments with a larger number of 
particles per cell (up to 32 per species), obtaining es- 
sentially the same results. Such a background of cold 
particles is initialized also within the current sheets, but 
with the addition of a hot component whose density pro- 
file is Uh = n/jo/cosh (. Here nho = rjUco is the den- 
sity of hot particles at the center of the current sheet. 
The temperature of this hot component is determined 
by pressure balance between the gaseous part (inside the 
current sheet) and the magnetic part (in the cold wind), 
which yields a thermal spread Qh = kTh/mc^ = o I2r\ 
for a 3D relativistic Maxwellian. Finally, the varia- 
tion of magnetic field across a current sheet needs to 
be sustained by a current flowing along z. In the wind 
frame, electrons and positrons in the current sheet should 
be drifting in opposite directions with the same speed 
^h = y/^ /(nio) (c/wp)/A. In summary, to initialize 
the distribution of hot electrons (and positrons) in the 
frame of the simulations, we need to make two succes- 
sive Lorentz boosts: first, from their proper frame to 
the wind frame, where electrons and positrons counter- 
stream along z with velocity ±/3;i; then, from the wind 
frame to the simulation frame, where the wind propa- 
gates along —X with Lorentz factor 70 |j 

In our study, we adopt 70 = 15 as our fiducial bulk 
Lorentz factor. As discussed in ^ we actually explore a 
wide range of Lorentz factors, from 70 = 3 to 70 = 375, 
and find basically the same results, modulo an overall 
shift in the energy scale. We vary the stripe wavelength 
by almost two orders of magnitude, from A = 20 c/wp to 
A = 1280 c/wp, and we discuss its effects on the structure 
of the shoc k and the process of shock-driven reconnec- 
tion in j ]4.1[ The dependence on the wind magnetization, 
though always in the regime of highly magnetized flows 
(cr 3> 1), is presented in ^4.2 where we vary a from 10 
to 100. Different values 01 th e stripe-averaged magnetic 
field are investigated in j ]4.3[ where we vary the param- 
eter a = 2{By)x/iBo + \{By)x\) between 0.0 and 0.95. 
For pulsar winds, large values of \a\ would be expected 
at high latitudes above the midplane, whereas a = 
along the equator (see upper panel in Fig. [I]) . 

The structure of the current sheet is completely de- 
termined once we fix the value of the density excess 
V = nho/uco, and of the sheet half-thickness A = X6/2tt. 
We choose rj = 3, motivated by analy tical studies of 
pulsar winds (Lyubarsky & Kirk 2001). However, we 
find that our results do not appreciably depend on the 
value of r] (we have also tried with 77 = 1.5 and 77 = 6), 
provided that the contribution of current sheets to the 
upstream particle flux stays negligible (i.e., rjS <?C 1). 
Regarding the sheet half-thickness, we typically employ 
A = c/ujp (for each choice of A, this determines the value 
of 5 = 27r A/A) , but we have verified that our results are 
not very sensitive to this parameter (we have tried up to 

^ In the range of parameters we explore, we find that /3;i <g; 1, 
so the first Lorentz boost can be safely neglected. 



A — Sc/wp, for A — 320c/tLip). The choice for the sheet 
thickness is guided, on the one hand, by the requirement 
that the contribution of current sheets to the upstream 
plasma energy flux should be small (i.e., (7(5 <C 1). On 
the other hand, we enforce A > O.bc/ojp, to ensure that 
the incoming current sheets are stable to the process of 
undriven reconnection. With this choice, we are confi- 
dent that the trigger for magnetic dissipation will come 
entirely from the shock-driven compression of the flow. 

We choose the relativistic electron (or positron) skin 
depth in the cold wind c/ujp such that the smallest scale 
in the system, which for ct > 1 is the relativistic Larmor 
radius r^ = (c/wp) Ax/ct, is resolved with at least a few 
computational cellsP] We resolve c/ujp with 7.5 compu- 
tational cells for a = 10, 10 cells for a = 20, 13 cells for 
a = 30, 17 cells for a = 50, and 24 cells for a = 100. For 
each value of the magnetization, we have verified that 
our results do not substantially change when doubling 
the resolution (for cr = 10, we have tried with up to 
c/ujp — 30 cells). Our computational domain is typically 
^ 400c/a;p wide (along the y direction). In particular, 
we choose a transverse size of 3072 cells for a — 10, of 
4096 cells for a = 20, of 5120 cells for a = 30, of 7168 
cells for a = 50, and of 9216 cells for a — 100. A large 
box is of paramount importance for the consistency of 
our findings. As we show in Appendix]^ the results of 
our simulations converge only when the computational 
box is larger than the half-wavelength of the stripes. It 
follows that ID sim ulatio ns, as the ones presented by 
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2007 1, cannot capture all the rele- 
vant physics, especially for long wavelengths. 

The results presented in the following sections refer to 
shocks that have already reached a steady state. For 
the typical parameters we explore, this happens before 
ojpt = 3000, but in some cases (especially for long stripe 
wavelengths), we need to evolve our simulations up to 
longer times. In Appendix [B} we follow the entire time 
evolution of a representative shock, from the earliest 
stages until it approaches a steady state. 

3. SHOCK STRUCTURE AND PARTICLE ACCELERATION 

In this section, we present the properties of a shock 
that propagates in a relativistic electron-positron flow 
with alternating fields (a "striped wind"). The case of 
uniform pre-s hock fields has been investigated by Sironi 
fc Spitkovsk7| (|2009) and |Sironi fc Spitkovsky] ( |2011j ) via 
multi-dimensional PIC simulations (tor electron-positron 
and electron-ion plasmas, respectively). They studied 
how the magnetic field strength and obliquity can affect 
the efficiency of the Fermi process, where particles gain 
energy by b ouncing back and forth across the shock (e.g.. 



Blandford & Eichler 1987). They found that in perpen- 



dicular shocks i.e., with held orthogonal to the shock nor- 
mal, the Fermi mechanism can operate only if the flow 
is weakly magnetized. For magnetizations a > 10~^, the 
field is so strong that charged particles are forced to slide 
along the field lines, whose orientation prohibits repeated 
crossings of the shock, thus inhibiting the Fermi process. 
In a medium with alternating fields, we find that the 
efficiency of Fermi acceleration is controlled by the mag- 

^ The skin depth and Larmor radius for the hot particles in the 
current sheet are larger than their counterparts in the cold wind 
by a factor of ^Jaj2rf- and cr/2ri, respectively. 



nitude of the stripe-averaged field {By)\. Striped flows 
with {By)\/ Bq <^ 1 are equivalent to weakly magnetized 
shocks in uniform fields, with respect to particle acceler- 
ation via the Fermi process. In this section, we choose a 
value for the stripe-averaged field {By)x/ Bq ~ 0.05 (cor- 
responding to a = 0.1) that is large enough to inhibit 
the Fermi process. This allows us to focus solely on the 
process of shock-driven reconnection, and to isolate its 
effects on the shock structure and the particle spectrum. 
In pulsar winds, the special case {By)\ = should be 
realized only along the equatorial plane of the wind. By 
choosing a = 0.1, we address the generic properties of 
the flow above (or below) the midplane. 

In this section, we adopt 70 = 15 and cr = 10 as our 
flducial values for the upstream bulk Lorentz factor and 
magnetization. To clarify the physics of shock-driven 
reconnection, we employ relatively large values for the 
stripe wavelength (A — 320c/wp or A = 640c/a;p), and 
we choose a = 0.1, corresponding to a stripe-averaged 
fleld {By)\/Bo ~ 0.05. The dependence of the shock 
properties on 70, c. A, and a is investigated in El There, 
we also explore the special case a — 0, where shock- 
driven reconnection coexists with Fermi acceleration. 

3.1. Shock Structure 

The typical structure of a relativistic shock propagat- 
ing in a striped wind is shown in Figs. [2] and [3j which 
refer to cOpt = 3750, after the flow structure has reached 
a steady state. The longitudinal profile of the transition 
region, along the direction x of shock propagation, is pre- 
sented in Fig.[2J whereas Fig.[3]shows the flow properties 
in the xy plane of the simulation, zooming in on a region 
around the shock (as delimited by the vertical dashed 
red lines in Fig. [2ja)). 

The incoming now {x > 2600c/ajp) carries an alternat- 
ing magnetic field of wavelength A = 640 c/wp (Fig. [2|c); 
black line for By, red for Ez). The field reverses its polar- 
ity across current sheets filled with hot overdense plasma 
(see the spikes in density and average particle energy of 
Fig. l2ja) and (d), respectively), where the magnetic en- 
ergy IS much lower than in the cold wind (see the troughs 
in the magnetic energy fraction eg — B'^ /STT^QmncoC^ , 
black line of Fig. ^h)). For x > 2600 c/wp, the only 
perturbation to the incoming flow comes from the so- 
called electromagnetic precursor, a coherent train of 
transverse electromagnetic waves propagating into the 
upstream at the speed of light (see the short-scale rip- 
ples at a; ~ 3850 c/wp in Fig. WCb) and (c)). The elec- 
tromagnetic precursor wave is generated in the initial 
stages of shock evolution, before shock-driven reconnec- 
tion appreciably changes the structure of the shock. At 
such early times, the flow resembles a magnetized shock 
propagating in uniform flelds (see Appendix pi), which 
is well kno wn to radiate electromagnetic waves into the 
upstream ([Gallant et al.||1992l IHoshinol 120081 ISironi fc 
2009||2011[ ). 



Spitkovsky 



i'he steady-state structure of the flow in Fig. [2] shows 
the presence of two shocks. The main shock (wnich we 
shall call the "hydrodynamic shock" from now on, for 
reasons that will become clear below) corresponds to 
the jump in density occurring at a; ~ 1000 c/wp. At 
X ~ 2600 c/wp, well ahead of the hydrodynamic shock, 
the incoming flow crosses a fast magnetohydrodynamic 
(MHD) shock, whose location is highlighted in Fig. [2] by 
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Fig. 2. — Internal structure of the flow at Upt = 3750, for stripe 
wavelength A = 640 c/ujp and magnetization cr = 10. The stripe- 
averaged field is {By)x/Bo ~ 0.05, corresponding to a = 0.1. The 
hydrodynamic shock is located at x ~ 1000c/i.ijp. In all panels, 
the location of the fast MHD shock (x ~ 2600c/a;p) is indicated 
with a vertical dotted blue line. As a function of the longitudinal 
coordinate x, the following quantities are plotted: (a) y-averaged 
particle number density n, in units of the upstream value; (b) y- 
averaged magnetic energy fraction eg = B^ /Sir-fomncoc^ (black 
line), and j/-averaged magnetic energy per particle eg/n (red line); 
(c) j/-averaged magnetic field By (black line) and electric field E^ 
(red line), normalized to the upstream magnetic field Bg; (d) mean 
kinetic energy per particle, in units of the bulk energy at injection; 
(e) X — 7/9x positron phase space; (f ) x — 7/3z positron phase space. 



a vertical dotted blue line. At the fast shock, the wind is 
decelerated and compressed, as shown by the modest in- 
crease in density (Fig. [2ja)), magnetic energy (black line 
in Fig.l^b)), and transverse magnetic field By (black line 
in Fig.Wc)). The deceleration of the flow appears most 
clearly m panel (c), as a mismatch between the motional 
electric field E^ and the magnetic field By behind the 
fast shock (red and black line, respectively). Irrespective 
of the initial bulk Lorentz factor of the flow, the plasma 
drift velocity (= —Ez/Byx) downstream from the fast 
shock is only moderately relativistic, since E^/By ~ 0.8 
there (Fig. l2jc)). The flow deceleration at the fast shock 
is accompanied by a partial isotropization of the parti- 
cle distribution, mostly in the xz plane orthogonal to the 
field. Overall, the flow kinetic energy decreases across the 
fast shock from 70 = 15 to (7) ~ 10 (Fig.[2]jd)), to com- 
pensate for the growth in magnetic energy (Fig. ^b)). 

The fast MHD shock has a dramatic impact on the 
structure of the striped wind. Across the fast shock, 
the stripe wavelength shrinks from A = 640 c/wp to 
A ~ SSOc/wp (see Fig. [2lc)). Most importantly, the 



compression induced by the passage of the fast shock 
through a given current sheet triggers the onset of mag- 
netic reconnection. Small-scale islands develop inside the 
current sheet as a result of the tearing-mode instabil- 
ity. The growth and evolution of reconnection islands is 
clearly shown in the 2D plots of density and magnetic en- 
ergy of Fig.jSJ^a) and (b), respectively. The current sheet 
at a: ~ 250(jc/wp, just downstream from the fast shock 
(located at a; ~ 2600c/ajp), appears partitioned into a 
number of overdense islands with strong magnetic fields 
(red point-like regions in Fig.lslb), at a; ~ 2500c/cL'p). A 
magnetic X-point exists in between each pair of neighbor- 
ing islands, where field lines of opposite polarity break 
and reconnect. As the flow recedes from the fast shock, 
the small-scale islands created inside each current sheet 
coalesce to form bigger islands. This temporal evolu- 
tion corresponds, in the snapshot taken at Wpt = 3750 in 
Fig. [SJ to a pattern of bigger (and so, fewer) islands for 
current sheets that are farther behind the fast shock. 

As reconnection proceeds with distance downstream 
from the fast shock, more and more magnetic energy is 
released and transferred to particles. In fact, the spikes 
in average particle energy of Fig. ^d) become stronger 
for current sheets farther from the JVIHD shock (see also 
the longitudinal phase space of positrons in Fig.[2Fe), for 
1000 c/wp "^ X < 2600 c/wp). The transfer of magnetic 
energy to the particles continues up to the point where 
reconnection islands fill the entire region between neigh- 
boring current sheets. Now, the striped structure of the 
flow is erased, and a hydrodynamic shock forms (located 
at X ~ 1000 c/wp in Figs.[2]and|3|. 

Behind the hydrodynamic shock, the flow comes to 
rest (in the wall frame of our simulations) , and the par- 
ticle distribution is isotropic in three dimensions (com- 
pare longitudinal and transverse phase spaces in Fig. ^e) 
and (f), respectively). The post-shock number density 
approaches rid — 4nu (Fig. ^a), where nu is the den- 
sity ahead of the fast shock), and the shock velocity is 
Psh — 1/3 (see Appendix [B|, in agreement with the jump 
conditions of a relativistic unmagnetized 3D plasma. In 
fact, most of the energy per particle downstream from the 
hydrodynamic shock is in kinetic form (as opposed to the 
dominant electromagnetic component of the incoming 
striped flow). This becomes clear when comparing the 
profile of the mean kinetic energy per particle (Fig.^fd)) 
with the ID plot of the mean magnetic energy per par- 
ticle (red line in Fig. ^b))rl Behind the hydrodynamic 
shock, the average kinetic energy per particle increases 
from (7)770 — 1 up to (7)770 ~cr-|-l~ 11, as expected 
in the case of full dissipation of magnetic fields. Corre- 
spondingly, the mean magnetic energy per particle drops 
to zero (red line in Fig. l2[b)). As best seen in experi- 
ments with a smaller stripe wavelength, in the case a 7^ 
a net magnetic field survives downs t ream from the shock, 
as predicted by Petri fc Lyubarsky (2007), and shown by 
the ID profile of By in i''ig.[2lc) (black line). This resid- 
ual magnetic field results from shock-compression of the 
stripe-averaged upstream field (we remind that a = 0.1 
here, corresponding to {By)\/Bo ~ 0.05). 

■* We point out that in the incoming flow electric and magnetic 
fields equally contribute to the energy balance, whereas the red line 
in Fig. pVo) only includes magnetic fields. However, downstream 
from thehydro dynamic shock, the flow is at rest (in the simulation 
frame), and electric fields vanish (see red line in Fig. [2[c)). 
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Fig. 3. — Internal structure of the flow at cjpt = 3750, for A = 640c/a;p, a = 10, and a = 0.1, zooming in on a region around the shock, 
as dehmited by the vertical dashed red lines in Fig. |2la). The hydrodynainic shock is located at x ~ lOOOc/cjp, and the fast MHD shock 
at a; ~ 2600c/ajp. The following quantities are plotted: (a) 2D plot in the simulation plane of the particle number density, in units of the 
upstream value, with contours showing the magnetic field lines; (b) 2D plot of the magnetic energy fraction eb = -B^/87r7omncOC^; (c) 
2D plot of the difference eg — e^, where the electric energy fraction is e^ = E'^ /Sn-fomncoc^', (d) 2D plot of the mean kinetic energy per 
particle, in units of the bulk energy at injection; (e)-(f) particle energy spectra, at different locations across the flow, respectively outside 
(panel (e)) or inside (panel (f)) of current sheets. The color of each spectrum matches the color of the corresponding arrow at the bottom 
of panel (d), showing where the spectrum is computed. The dotted line in panel (e) is a Maxwellian with the same average energy as the 
downstream particles; the dashed lines in panels (e) and (f) indicate a power-law distribution with slope p = 1.4. 



We point out that the picture described above rep- 
resents accurately the long-term behavior of the shock. 
In particular, the distance between the fast and the hy- 
drodynamic shock (~ 1600 c/ujp in Figs. [2] and ^ stays 
approximately constant in time (see Appendix Bl). The 
hydrodynamic shock moves with f3sh ~ 1/3, and the fast 
MHD shock remains far enough ahead of the main shock 
such that to guarantee that reconnection islands will fill 
the space between neighboring current sheets, by the 
time the flow enters the hydrodynamic shock. For a fixed 
reconnection rate, this implies that the distance between 



the fast and the hydrodynamic shock will scale linearly 
with the stripe wavelength (as will the size of reconnec- 
tion islands just ahead of the main shock), a prediction 
that we have directly verified in our simulations. 

Finally, we refer to Appendix [K\ for a detailed discus- 
sion of the differences between our 2D results and the 



ID simulations presented by Petri & Lyubarsky (2007). 
As apparent in Fig. [3J the process ot island coales- 
cence, which is essential for the formation of the hy- 
drodynamic shock, can be captured correctly only with 
multi-dimensional simulations. Indeed, for the parame- 



ters employed in this section, the ID model by [Petri "^ 
Lyubarsky (20071 would predict negligible field dissipa- 
tion, in sharp contrast with our findings. This clearly 
emphasizes the importance of multi-dimensional physics 
for our understanding of shock-driven reconnection. 

3.2. Particle Spectrum and Acceleration 

The particle energy spectrum at different locations 
through the flow, as marked by arrows at the bottom 
of Fig. Isjd), is shown in Fig. Isle) and (f). Panel (e) 
shows the particle distribution in the cold wind, whereas 
panel (f ) focuses on the hot plasma within current sheets. 
Downstream from the hydrodynamic shock, the striped 
structure is completely erased, and no difference persists 
between current sheets and cold wind. Here, the particle 
spectrum (black line in both panels) is in the form of a 
flat power-law tail (with spectral index p ~ 1.4, dashed 
line) extending from 7ijiin ~ 30 to 7max — 500, where 
it cuts off exponentially. Since 1 < p < 2, most of the 
particles are found at 7 ~ 7min, but most of the energy 
is contributed by particles at 7 ~ 7max- For comparison, 
a 3D Maxwellian with the same average Lorentz factor 
(7) = 7o(f + 1) — 165 is plotted in panel (e) as a dotted 
line, to show that the actual particle spectrum is much 
broader than a thermal distributionr] By following the 
flow from the fast to the hydrodynamic shock, we now 
clarify how such a particle spectrum is generated. 

Most of the incoming plasma is contained in the cold 
striped wind exterior to current sheets, and its spectrum 
is initially a cold Maxwellian drifting with bulk Lorentz 
factor 7o = 15 (light blue line in Fig. [3]je)). After cross- 
ing the fast shock, the flow decelerates and it partly 
isotropizes, which explains why the particle spectrum is 
now broader and it peaks at lower energies (dark blue line 
in Fig. l3|e) , but see also the 2D plot of average particle 
energy m Fig. [3ld), and its ID projection in Fig. [21(d)). 

The evolution of the particle population within current 
sheets is much more dramatic, as shown in Fig. p[f)r| 
Ahead of the fast shock, the spectrum at x ~ 2850c/ajp 
(yellow line in panel (f)) shows, in addition to the cold 
plasma in the striped wind, the hot particles in the cur- 
rent sheet. Their distribution is initialized as a drift- 
ing 3D Maxwellian with thermal spread Qh = o'/2?7, so 
that in the simulation frame the spectrum should peak 
at (7) = 37oO/i ~ 75, as observed (we remind that 
cr = 10 and 77 = 3). Behind the fast shock (red line 
for X ~ 2200 c/wp, brown line for x ~ 1350 c/wp), the 
particle spectrum consists of two components: a low- 
energy peak, reminiscent of the particle distribution far 
from current sheets (compare with the dark blue line in 
Fig. 3[e)); and a high-energy component, that grows in 
number and extends to higher and higher energies, as the 
flow propagates toward the hydrodynamic shock. This 
component is populated by particles that were initially 
outside the current sheet, and have entered the sheet 
in the course of the reconnection process. They mostly 
end up in the reconnection islands of Fig. [31 that act as 
reservoirs of particles (Fig.lsFa)) and particle energy (the 

^ We remark that the particle spectrum behind the hydrody- 
namic shock does not significantly evolve with time, after the shock 
has reached a steady state (see Appendix ml. 

® We remind that the spatial pattern shown in Fig. [sl as a func- 
tion of distance behind the MHD shock, corresponds to the time 
history of a given fluid element after its crossing of the fast shock. 



high-energy regions of panel (d), colored in red, match 
very well the overdense islands of panel (a)). In fact, the 
main contributions to the high-energy end of the spec- 
trum at X ~ 1350 c/wp (brown line in Fig. ^f)) come 
from the two major islands of that current sheet. 

As magnetic reconnection proceeds with distance be- 
hind the fast shock, particles are accelerated to higher 
and higher energies, as we discuss below in more de- 
tail. This continues until reconnection islands grow as 
big as the distance between neighboring sheets, at which 
point the current sheets occupy the entire flow (and the 
hydrodynamic shock forms). This explains the smooth 
transition from the spectrum within current sheets to the 
distribution downstream of the hydrodynamic shock (yel- 
low through black lines in Fig.l3[f)). Also, since most of 
the particles that end up with high energies were initially 
part of the cold wind, the resulting spectrum behind the 
hydrodynamic shock is almost independent of the initial 
setup of the current sheet (most notably, sheet thickness 
A and overdensity 77). 

We now investigate in more detail the physics of parti- 
cle acceleration within current sheets. In Fig.|4J we follow 
the trajectories of a representative sample of positrons, 
extracted from the simulation of a striped wind with 
A = 320c/a;p (but our conclusions hold for all wave- 
lengths). All the selected positrons are initially in the 
cold wind (i.e., far from current sheets), at roughly the 
same x-location (within a range of 5c/u;p). At first they 
propagate toward the fast shock, which they cross at 
iOpt ~ 700, as shown by the spreading of their energies in 
Fig. Wla,). They encounter the closest current sheet (in 
this case, the one that was to the right of their initial 
position) at uipt ~ 1008, as shown in Fig. W[h) together 
with the corresponding 2D plot of particle number den- 
sity. Starting from this time, their histories can diverge 
significantly, depending on their y-location at the mo- 
ment of interaction with the current sheet. The particles 
that will end up with relatively low energies (7 < 370, 
black lines in panel (a)) are found at ujpt = 1008 exclu- 
sively around magnetic islands (their locations are indi- 
cated with filled circles in panel (b)). At later times, they 
will get trapped on the outskirts of the growing islands, 
without appreciable changes in energy. 

On the contrary, the positrons that will eventually 
reach high energies (7 > 670, red and blue lines in 
Fig. Ba)) are concentrated at Wpi — 1008 in the vicin- 
ity ofmagnetic X-points (the particle locations are shown 
with open circles in Fig.l4Fb)). It is at X-points that mag- 
netic field lines of opposite polarity are converging and 
tearing (white contours in Fig. WCh) show the magnetic 
field lines), and magnetic energy is transferred to parti- 
cles. Starting from ujpt — 1008 (vertical dashed line in 
panel (a)), the energies of these positrons grow simulta- 
neously and explosively (red and blue lines in Fig. |4l^a)), 
due t o the reconnection electric field parallel to tne X- 
line ([Zenitani fc Ho sh"ino"2 001| paroschek et al.[[2004[ 
Lyubarsky & Liverts 2008 ) . We find that X-point ac- 
ceieration represents the first stage of energization for 
most of the particles that will end up with high ener- 
gies. Particles pre-accelerat ed at X-points continue to 
gain energy by compression ( Lyubarsky &: Liverts[[2008"' 



Giannios 2010), as they are advected from the X-pomt 
into the closest island. No major systematic changes of 
energy are observed after the particles get trapped within 
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Fig. 4. — Energy evolution of a sample of representative positrons 
interacting with a current sheet, in a striped flow with A = 
320c/a;p, a = 10, and a = 0.1. All the selected positrons have been 
initialized as part of the cold wind, at roughly the same x-location. 
Panel (a): Energy histories of the selected positrons, with lines of 
different color depending on their final energy (black if 7 < 870, 
blue or red if 7 > 670) and their j/-location at ujpt = 1008 (blue if 
y < lOOc/tjp, red if y > lOOc/ojp). Panel (b): Magnetic field lines 
(white contours) superimposed on the 2D plot of particle number 
density, at time LOpt = 1008 (marked as a vertical dashed line in 
panel (a)); the locations of the selected positrons are shown as yel- 
low filled or open circles, depending on their final energy (filled 
circles for particles with 7 < 870, open circles for 7 > 670). 



magnetic islands. Different X-points may provide differ- 
ent energy gains, as shown by the red (respectively, blue) 
curves in panel (a), which refer to high-energy positrons 
that get accelerated at the X-point above (respectively, 
below) the small magnetic island in the middle of the 
current sheet. In summary, the energy gain of a given 
particle is determined on the one hand by the strength 
of the reconnection electric field at the X-point, and on 
the other hand by the time available for acceleration, as 
the particle flows from the X-point into the closest island. 
In retrospect, the importance of X-points for the gen- 
eration of nonthermal particles is not surprising. In 
magnetically-dominated flows, particles can untie from 
magnetic field lines and be injected into the acceleration 
process only in regions where the electric field is compa- 
rable or exceeds the magnetic field, i.e., green or blue re- 
gions in the current sheets of Fig. ISTc) (where we plot the 
difference es — ^E between magnetic and electric energy). 
Based on this criterion, X-points are natural sites of par- 
ticle energization, since the magnetic field vanishes there. 
As discussed above, the main agent of particle accelera- 
tion at X-points is the reconnection electric field, which 
for the magnetic geometry employed in this work is ori- 



ented along z (more precisely, along +z if the magnetic 
field By is negative to the left of the current sheet and 
positive to the right; along — z in the opposite case). The 
sign of the reconnection electric field leaves an imprint 
on the particle z-momentum. As shown by the positron 
transverse phase space in Fig. [2jf), current sheets where 
the reconnection field is expected to be along +z are 
mostly populated by high-energy positrons with 7/3z > 
(and electrons with j(3z < 0), and the opposite holds 
where the reconnection field is along —z. 

Accelerated by the reconnection electric field, high- 
energy particles stream at first along the z direction, 
and then, by effect of the reconnected magnetic field 
Bx, they are advected along y into the closest island. 
In our 2D simulations, where the z dimension is not re- 
solved, the relevant timescale for acceleration is the time 
to drift from an X-point into the closest magnetic island. 
As the fiow propagates downstream from the fast shock, 
magnetic islands grow and merge (see Fig. [sla)), so that 
the typical distance between an X-point and the closest 
island increases, allowing particles to be accelerated to 
higher energies. This explains why the upper cutoff in 
the particle spectra of Fig. IsFf ) shifts to higher energies 
for current sheets farther from the fast shockQ 

Since our 2D simulations cannot resolve potential per- 
turbations to the particle motion along z, the validity of 
our conclusions may be questioned, if the structure of X- 
points in a realistic 3D scenario is much different from the 
2D picture analyzed here. In particular, in 3D we expect 
the current sheet to be folded, by effect of the so-calle d 
drift-kink instability (Zenitani & Hoshino 2005a 2007). 
The characteristic wavelength ot the kink mode in the xz 
plane will introduce a different length scale, that could 
potentially compete with the distance between X-points 
and islands (in the xy plane) in constraining the max- 
imum energy of accelerated particles. In fj5] we address 
this important point, and show that the latter constraint 
is more restrictive than the former, at least for the pa- 
rameters explored in this work. It follows that our 2D 
simulations with in-plane fields can capture very well the 
main 3D properties of the shock, and in particular the 
physics of particle acceleration. 

Although we find that the reconnection electric field 
at X-points is the main agent of particle acceleration, 
we have verified, by tracing the orbits of a large num- 
ber of particles extracted directly from the simulation, 
that a variety of other mechanisms can play a role in 
particle energization. During the coalescence of mag- 
netic islands, particles may be accelerated by the anti- 
reconnection electric field at the X-point located in be- 
tween the two m erging islands (jJaroschek et al. 2004 
Oka et al.| |2010|) . Alternatively, particle energization 



may be due to a kind of first-order Fermi process within 
magnetic islands, when particles reflect fr om the ends 
of con tracting islands following coalescence (Drake et al. 
2006|. Also, the second-order Fermi process can operate 



between randomly-moving islands just downstream from 
the hydrodynamic shock. The two latter mechanisms 
change primarily the components of particle momentum 

^ Another factor that drives the trend in the upper spectral cutoff 
of Fig. ISlf), though less important than the timescale argument 
mentioned above, is the increase in the strength of the reconnection 
field as the flow is compressed toward the hydrodynamic shock. 



along the plane where the magnetic field lies (the xy 
plane, in our simulations). On the contrary, for recon- 
nection and anti-reconnection at X-points, it is the mo- 
mentum component along the reconnection electric field 
(along z, in our case) that receives the greatest boost. 

The energization mechanisms discussed so far are pe- 
culiar to the process of magnetic reconnection, being 
seeded by the field dissipation itself, or by the resulting 
fluid motions. In shock-dominated systems, two more 
acceleration processes can be relevant, whose free en- 
ergy comes from the converging motion of the pre-shock 
and post-shock flows. In first-order Fermi acceleration 
(or diffusive shock acceleration, DSA), particles diffuse 
back and forth across the shock front and gain energy 
by scattering from magnetic tur bulence embedded in the 
converging flows (e.g., [B landford fc Ostriker||1978 

— ;E 



1978 Blandford & EicI 



Bell 

ier||1987p. In magnetized flows, 



particles may also gain energy directly from the large- 
scale motional electric field while they gyrate around the 
shock front , a process known as shock-drift acceleration 
(SPA; e.g.,|Chen fc ArmstronK||1975[|Webb et al. 



1983 



Begelman fc KirE[l99l)[ ). As we describe in the toUow 
mg section, the JSljA process is most efficient in shocks 
where {By)x/BQ > 0.015, whereas the DSA mechanism 
can only operate if the stripe-averaged field is small, 
{By)x/B,, < 0.015. 

4. DEPENDENCE ON THE WIND PROPERTIES 

In this section, we investigate how the physical condi- 
tions in the wind can affect the structure of the shock and 
the physics of magnetic reconnection, with particular fo- 
cus on the particle spectrum downstream from the hydro- 
dynamic shock. We will explore th e de pendence of our 
results on the s trip e wavelength (in ^ 4.1 ), the wind mag- 
netizatio n (in [4.2), the strength of the stripe-ave rage d 
field (in { 4.3 ), and the wind b ulk Lorentz factor (in J4.4). 

The ID model presented by Petri & Lyubarsky (20071 
predicts that magnetic dissipation should not appre- 
ciably affect the dynamics of the flow, provided that 



Ac/^„/<^- 8Ci/5. Here, A, 



./^ — X/{c/ujp) is the stripe 
wavelength in units of the relativistic skin depth in the 
cold wind, and fi = 6 — 10 is the ratio of downstream 
sheet thickness to downstream Larmor radius (its value 
is calibrated with ID PIC simulations). In this case, the 
striped structure of the wind will be preserved down- 
stream, just compressed, and the flow will satisfy the 
MHD jump conditions for a strongly magnetized fluid. 
On the contrary, full dissipation of the alternating flelds 
should occur if Xc/uj /V^ ^ 4^i/5, when the thickness 
of downstream current sheets becomes comparable to 
the distance between ne ighboring sheets. In this case, 
Petri & Lyubarsky (2007) flnd that the post-shock flow 
is purely hydrodynamicai, since all the magnetic held en- 
ergy has been transferred to the particles. 

As we discuss below, the results of our 2D and 3D 
simulations suggest a different picture. We find that 
complete dissipation of the alternating fields by magnetic 
reconnection is a common by-product of the shoc k evo- 
lution, eve n in the cases when the ID model of |Petri| 
& Lyubarsky (2007) predicts that the stripes should be 
preserved downstream from the shock. As suggested by 
Fig. |3l multi-dimensional simulations are of key impor- 
tance to correctly capture the evolution of the tearing- 
mode instability, which is ultimately responsible for field 



dissipation and particle acceleration. 

Even though complete field dissipation (and transfer 
of field energy to the particles) occurs irrespective of the 
wind properties, the shape of the downstream particle 
spectrum is sensitive to the stripe wavelength and the 
wind magnetization. In most cases, the particle spec- 
trum can be described as a power law in energy, with a 
flat slope of index 1 < p < 2. The slope is moderately 
sensitive to the pre-shock parameters, but it is the width 
of the spectrum that depends most dramatically on the 
properties of the wind. If the spectrum can be modeled 
as a power law of index p between 7mi„ and 7niaxi we 
can derive a relation between the lower and upper cut- 
offs by requiring that the mean particle Lorentz factor 
be (7) = 7o(tT -I- 1), as expected for full dissipation of 
magnetic energy. We obtain 
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For 1 < p < 2 and fixed a, we see that this expres- 
sion includes both the case of narrow distributions with 
7niin ~ 7max — 700": e^h^ the casc of wide spectra with 
7o and 7 



ri/(2-p) 



4.1 



7min - 7o ana 7max ^^70^""' ^'. As we show in 
the stripe wavelength is the main factor, at fixed wine 
magnetization, that controls the transition from narrow 
Maxwellian-like spec tra to broad power-law tails. 

In Figs. 5][7 9] and 12 we show the dependence of the 
downstream particle spectrum on the conditions in the 
wind. In most cases, the spectrum is measured between 
700 c/ujp and 300 c/wp downstream from the hydrody- 
namic shock|^ In the subpanel of each figure, we es- 
timate the mean particle kinetic energy (7)/7o (black 
line) and the width of the particle spectrum, measured 
as the ratio 7max/7min of upper to lower energy cutoff 
(red line). More precisely, 7mi,i will be defined as the 
Lorentz factor where most of the particles lie (i.e., where 
jdN/d'y peaks), whereas 7niax will correspond to the 
location where most of the energy resides (i.e., where 
j'^dN/dj is at maximum). 

4.1. Dependence on the Stripe Wavelength 

Fig. [5] shows how the downstream particle spectrum 
changes for different stripe wavelengths A, keeping fixed 
the magnetization a = 10, the stripe-averaged field 
{By)x/BQ ~ 0.05 (corresponding to a = 0.1), and the 
bulk Lorentz factor 70 = 15. We explore a wide range of 
wavelengths, from A = 20c/a;p up to A = 1280 c/wp. 

Regardless of the value of A, the s truct ure of the flow 
resembles the picture presented in |3.1| A fast MHD 
shock propagates into the striped wind, compressing 
the incoming current sheets and triggering magnetic re- 
connection. Due to reconnection, the striped structure 
is eventually erased, and a hydrodynamic shock forms. 
Downstream from the shock, the flow is nearly unmag- 
netized, apart from the field component associated with 
the shock-compression of the stripe- averaged field {By)\. 
It follows that the downstream mean kinetic energy per 
particle should be close to the value (7)770 = (J + 1 ex- 
pected in the case of complete field dissipation. The value 

* The width of the slab where we compute the spectrum is cali- 
brated such that to average out all potential variations along x, and 
to provide an estimate of the characteristic downstream spectrum. 
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Fig. 5. — Downstream particle spectrum at LOpt = 3000 for dif- 
ferent stripe wavelengths A, in a flow with cr = 10 and a = 0.1. We 
vary the stripe wavelength from A = 20c/ujp up to A = 1280c/ajp. 
The dotted red line is a Maxwellian with the same average en- 
ergy as the spectrum colored in red (which refers to A = 20c/ujp); 
the dashed black line indicates a power-law distribution with slope 
p = 1.5. In the subpanel, the black line shows the average Lorentz 
factor of downstream particles (7) as a function of the stripe wave- 
length A (axis on the left side). In the same subpanel, the red line 
presents the dependence of 7max/7min on ^ (axis on the right side). 
Here, 7miii is defined as the location where ■ydN/d'y peaks (i.e., 
where most of the particles reside), whereas 7max is the Lorentz 
factor where 'y-^dN/d-y peaks (i.e., where most of the energy lies). 



of (7) measured in our simulations (black line in the sub- 
panel of Fig.|5| is approximately constant with respect to 
A, and is consistent with such expectati onp] Clearly, this 
is in c ontrast with the ID model of Pet ri fc Lyubarsky] 
(2007), which would predict negligible field dissipation if 
'p (for ^1=8 and a — 



A ^T30 c/ujp (for ^1 = 8 and cr = 10). 

Even though the mean kinetic energy per particle does 
not appreciably vary with A, the shape of the spectrum 
does change, with a clear tendency for broader spectra 
at longer stripe wavelengths. As shown by the red line in 
the subpanel of Fig. pi the ratio 7max/7min increases with 
wavelength from 7max/7min — 1 (at A = 20c/ujp) up to 
7max/7min — 20 (at A = 1280 c/wp). Longer wavelengths 
yield smaller values of 7min (see main plot), which then 
imply larger values of 7max, given the relation in eq. (pi). 

This trend can be easily understood by considering the 
structure of t he fl ow for different stripe wavelengths. As 
described in §3.1[ the hydrodynamic shock is located at 
the point where reconnection islands grow large enough 
to fill the entire space in between neighboring current 
sheets. Since the distance between sheets is proportional 
to A, the same scaling should hold for the size of re- 
connection islands, just upstream of the hydrodynamic 
shockj^ Shorter wavelengths will then result in more 

® The trend of smaller values of ('y)/-fo with decreasing A (black 
line in the subpanel of Fig. pjl is due to the fact that the fractional 
contribution of current sheets (with fixed thickness 2 A ~ 2 c/ujp ) 
to the upstream flow is larger for smaller A. When averaged over 
one stripe wavelength, this results in a lower upstream energy per 
particle for smaller A (including both magnetic and kinetic contri- 
butions), which is then reflected in the post-shock value of {■y)/-fo- 

^" In the remaining of this section, by "reconnection islands" we 
will be referring only to the islands in the current sheet just ahead 
of the hydrodynamic shock, where the shape of the downstream 



numerous islands of smaller size, whereas fewer but big- 
ger islands will be present for longer A (see Fig. Isl . Since 
an X-point exists in between each pair of neighboring 
islands (belonging to the same current sheet), the num- 
ber of X-points per unit length (along the current sheet) 
will be larger for smaller A. For short wavelengths, most 
of the incoming particles will likely pass in the vicin- 
ity of one of the numerous X-points, as the flow crosses 
the hydrodynamic shock. The energy evolution of one 
particle will then be similar to that of any other parti- 
cle, with its Lorentz factor increasing from 7 ~ 70 up 
to 7 ~ 7o(cr + 1), as the particle gains energy from the 
annihilating fields. This explains why for short stripe 
wavelengths (A = 20c/a;p, red curve in Fig. I5) the par- 
ticle spectrum is so narrow, and similar to a Maxwellian 
distribution (red dotted curve). In summary, for short 
stripe wavelengths, all particles are equally close to an 
X-point, so they gain comparable amounts of energy. 

On the other hand, for long wavelengths, the energy 
evolution of diffe rent particles can be extremely diverse, 
as discussed in ^3.2| Particles that stay away from X- 



points are likely to remain cold, and retain the energy 
7 ~ 7o they started with. On the contrary, particles that 
pass through an X-point can gain a significant amount 
of energy. With increasing stripe wavelength, more par- 
ticles will belong to the former group, and fewer to the 
latter, since the number of X-points decreases. This ex- 
plains why the peaks in the spectra of Fig. [5J which track 
the downstream Lorentz factor of the "typical" particle, 
shift to lower energies with increasing stripe wavelength. 
For the same reason, the best-fitting power law becomes 
softer, with spectral slope increasing from p ~ 1.0 to 
p ~ 1.5 (indicated as a black dashed line in Fig. [sl. In 
response to a lower 7inin and a steeper slope, the up- 
per cutoff of the spectrum shifts to higher energies, as 
predicted by eq. ([2]). This is just a consequence of the 
longer time available for acce lerat ion, at larger stripe 
wavelengths. As discussed in j ]3.2[ particles injected at 
X-points are continuously accelerated by the reconnec- 
tion electric field, while they drift from the X-point into 
the closest island. Since the spacing between neighboring 
islands (in the same current sheet) increases with stripe 
wavelength, the injected particles will be able to stay in 
the acceleration region for longer times (and so, to reach 
higher energies), if A is la rger. 

As we discuss in ^4.2 the threshold between short 
and long wavelengths depends on the wind magnetiza- 
tion. Here, by "short A" we generically mean all cases in 
which the spectrum is narrow, similar to a Maxwellian. 
In contrast, the cases with "long A" have broad power- 
law spectra with slope 1 < p < 2. W e find that the con- 
dition proposed for full dissipation by [Petri fc Lyubarsky] 
(2007), namely that Xc/wp/V'^ ^ 4^i/5, works fairly 
well to discriminate between short and long wavelengths. 
For ^1=8, such criterion would predict Maxwellian-like 
spectra for A < 20c/i:jp (if a = 10), A < 45c/a;p (if 
a = 50), and A < 65c/a;p (if a = 100), which is indeed 
what we observe in our simulations. In retrospect, this is 
not surprising. In fact, the condition Xc/uj I \f^ ^ 4^i/5 
comes from requiring that the downstream sheet thick- 
ness be larger than the separation of neighboring sheets. 
In this limit, no distinction should persist between the 

particle spectrum is established, as demonstrated by Fig. pff). 
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plasma within (versus outside) the current sheets, i.e., all 
particles should share the same energy evolution. This 
obviously results in a Maxwellian distribution. 

Finally, we point out that the minor bump emerging 
at high energies in the spectra of Fig. [5] (most notably 
for wavelengths A < 160c/wp) is populated by parti- 
cles that are energized via the shock-drift mechanism 
(SDA) at the hydrodynamic shock. Such particles gain 
energy from the stripe-averaged motional electric field 
{Ez)\ — {By)\, while they gyrate around the shock. 
Since {Ez)\ > for our simulation setup, shock-drift 
accelerated positrons should preferentially have 7/3z > 0, 
whereas the opposite should hold for electrons. The few 
high-7/3z positrons seen at a; ~ llOOc/wp in panel (f) of 
Fig. [2] are indeed accelerated via the SDA mechanism. 
They gyrate around the shock a few times, before being 
advected downstream by the stripe-averaged magnetic 
field. Given the limited number of acceleration cycles, 
the SDA component in the spectra of Fig.[5]does not ex- 
tend in time to higher energies. So, SDA is not a promis- 
ing candidate to produce broad power-law tails. 

The efficiency of injection into the SDA process is 
higher for shorter stripe wavelengths, as Fig. [5] suggests. 
As seen from the pre-shock frame, a particle coming 
from the downstream side can be accelerated via SDA 
in the upstream stripe-averaged field only if it can sam- 
ple the full wavelength of the striped wind during its 
upstream residence time. Since particles returning up- 
stream are caught up by the shock after having trav- 
elled a distance r^u/lo (here, r^ u is the particle Lar- 
mor radius measured in the upstream frame) , we require 
7o(A/2) < rL,n/lo for efficient SDA. If 7 = ^2(7) - 67oo- 
is the Lorentz factor of particles at the high-energy end 
of the distribution (^2 = 1 — 10, depending on A), the 
previous condition can be rewritten as Xc/u /v^ ^2^2- 
For fixed magnetization, the criterion for efficient SDA is 
more easily satisfied at short stripe wavelengths, which 
explains why the normalization of the high-energy bump 
in Fig. [5] increases for smaller A. 

We remark that, apart from the factor ^2, which could 
in principle depend (though weakly) on A and a, the cri- 
terion for efficient SDA involves the same combination 
of wind magnetization and stripe wavelength (namely, 
Ac/w /v^) ^s the condition discussed at the beginning 
of this section, which regulates how closely the parti- 
cle spectrum resembles a Ma xwellian distribution (see 
also Petri & Lyubarsky 20071. Since Xc/^ is the wind 



A = 160 c/cjp a = 0.1 



cjpt=3000 



wavelength normalized to the skin depth of the pre-shock 
fluid, and since the average particle energy increases by 
a factor of ~ cr across the shock as a result of field dissi- 
pation, then Xc/uj /V^ is just the wind wavelength mea- 
sured in units of the post-shock plasma skin depth (apart 



from factors of order unity). In [4.2 we further com- 
ment on the importance on the parameter A,,/^ /V^ i^ 
regulating the physics of the shock and the shape of the 
downstream particle spectrum. 

4.2. Dependence on the Wind Magnetization 

In this section, we investigate the dependence of our re- 
sults on the magnetization of the wind. Fig.lHlshows how 
the downstream spectrum changes with a, keeping fixed 
the stripe wavelength A = IGOc/wp, the stripe-averaged 
field {By)\ ~ 0.05 (corresponding to a = 0.1), and the 
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Fig. 6. — Downstream particle spectrum at uj^t = 3000 for dif- 
ferent magnetizations c, in a flow with A = 160c/ajp and a = 0.1. 
We vary the wind magnetization from cr = 10 to cr = 100. In the 
subpanel, the black line shows the average Lorentz factor (7) of 
downstream particles as a function of cr (axis on the left side). The 
value of (7) is normalized to the total energy per particle (kinetic + 
electromagnetic) in the pre-shock flow, i.e., 700-. In the same sub- 
panel, the red line presents the dependence of the spectral width 
7max/7min On o" (axis on the right side). 



bulk Lorentz factor 70 = 15. We explore a wide range of 
magnetizations, from ct = 10 up to cr = 100. 

We find that efficient annihilation of magnetic fields by 
shock-driven reconnection, and transfer of magnetic en- 
ergy to the particles, occur irrespective of the wind mag- 
netization. As shown by the black line in the subpanel of 
Fig. [6j the average Lorentz factor of downstream parti- 
cles approaches in all cases the value (7) ~ 70CT expected 
for complete field dissipation. As the wind magnetization 
increases, the spectrum shifts to higher Lorentz factors 
(just because (7) oc a) and it changes in shape, with the 
part at low energies getting de-populated at the expense 
of the high-energy component. As a result, the particle 
spectrum, which could be fitted as a broad power law 
of index p ~ 1.2 for cr = 10 (blue curve), approaches a 
Maxwellian- like distribution for a = 100 (red curve). As 
a consequence, the ratio 7max/7min, a proxy for the width 
of the spectrum, becomes smaller for increasing a (red 
line in the subpanel of Fig. pi with axis on the right). 

This is the same trend observed when decreasing A 
at fixed a, as discussed in j ]4.1[ In the limit of very high 
magnetizations, for fixed wavelength (or very short A, for 
fixed cr), the downstream energy spectrum approaches a 



Maxwellian distribution, when A , 



■.lu 



o'^ 4^1/5 (here. 



^1 = 6 - 10). For A = 20c/a;p, 40c/wp, and SOc/wp, we 
have verified that the value of a above which the spec- 
trum resembles a Maxwellian is in good agreement with 
this criterion. For A — IGOc/wp, we would expect a 
Maxwellian distribution only for a > 300 (if ^1 = 8), 
beyond the range of magnetizations we explore. How- 
ever, the trend toward a single-component distribution 
(as a increases) is already clear within the limited range 
of magnetizations investigated in Fig. |6] 

The dominance of the high-energy component at large 
magnetizations can be understood in terms of the basic 
properties of the flow. In experiments with fixed a and 
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different A (see ^.ll, we have related the shape of the 



,1/3 



spectrum to the size of magnetic islands just ahead of 
the hydrodynamic shock, that scales with stripe wave- 
length. Here A is fixed, and the change in spectral 
shape with magnetization should instead be attributed 
to the different reconnection rate, namely the velocity 
/3rec at whi c h cold plasma flows into the X-points. |Zen-| 



^c /Vcr =160/50'^' a = 0.1 cjp(50/cr)^/H=300 



1/3. 



itani et al. ( 2009] ) have measured the reconnection rate 
in high-cr pair plasmas via 2D two-fluid MHD simula- 
tions of undriven reconnection, finding that the inflow 
velocity scales roughly as /^j-oc — 0.07-yCT in the range 
cr = 10— 100. We have confirmed this scaling in our sim- 
ulations, by measuring as a function of magnetization (at 
fixed A) the distance between the fast and the hydrody- 
namic shock. Everything else being fixed, this distance 
should be inversely proportional to the reconnection rate, 
and in the range of magnetizations explored in our study 
we find that the reconnection speed scales as jS-^^c oc \f^- 

As the reconnection rate increases with magnetization, 
a larger number of particles can be advected into a given 
X-point per unit time, which results in more particles ex- 
periencing the reconnection electric field. In other words, 
at larger a it is more likely for the "typical" particle to 
enter an X-point and gain energy from field dissipation. 
This explains why for high magnetizations the spectra 
in Fig. [6] are dominated by a single component of hot 
particles accelerated by the reconnection field. A similar 
argument explains why, for large values of cr, the spec- 
trum does not extend much beyond the average Lorentz 
factor (7). In response to the increase of inflow speed 
with a^ also the outflow velocity from a given X-point 
into the closest island increases with magnetization. For 
fixed A, which corresponds to a fixed spacing between 
X-points and magnetic islands, the time available for ac- 
celeration by the reconnection field will be shorter, for 
larger a. So, the maximum energy at which particles 
can be accelerated will be lower (relative to the average 
Lorentz factor (7)), for higher magnetizations. 

Finally, we notice that the minor bump emerging at 
high energies from the spectra of Fig. [61 which is pop- 
ulate d by particles accelerated via the sDA mechanism 
(see ^4.1), becomes more prominent for larger magne- 
tizations. This is in good agreement wit h th e criterion 




Fig. 7. — Downstream particle spectrum for different values of 



A and a, but such that the ratio A, 



,/,. 



A 



c/u. 



■ 160/V50 ~ 22.5 is 
fixed (as well as o = 0.1). Here, X^/^ is the wind wavelength 
in units of the skin depth in the pre-shock fluid, so Xc/ui /v^ i^ 
roughly the wavelength normalized to the skin depth of the post- 
shock fluid. The comparison is performed at the same time, in 
units of the relativistic plasma frequency of the post-shock flow (~ 
u)p/ y/a). Spectra are shifted along the x-axis by 50/(t to facilitate 
comparison with the case cr = 50. The black line in the subpanel 
shows the average downstream Lorentz factor, as in Fig. [6] 

the reference case cr = 50 (black curve). This compen- 
sates for the overall shift to higher Lorentz factors due 
to the increase in (7) ex cr, as expected for complete field 
dissipation!^ We confirm that the shape of the spec- 
trum, at thenigh-energy end, is controlled by the single 
parameter A c/ui I \f^- For a fixed value of this particular 
combination of A and a^ both the high-energy component 
of the spectrum and the bump of shock-drift accelerated 
particles are insensitive to variations in either A or ct. 

The shift of the low-energy part to higher Lorentz fac- 
tors, as a decreases, is just an artificial consec^uence of 
our re-scaling in the a;-axis of Fig. [Tj In reality, the low- 
energy cutoff stays fixed at some multiple of the pre- 
shock bulk Lorentz factor 70, irrespective of A or a, pro- 
< 2 ^2 discussed at the end of 5J4T| which die- vided A ^.j^^ j sjo is kept constant. With decreasing A, the 

number of X-points per unit length increases as oc A~^, 
which would make it easier for the "tyipical" particle to 
gain energy from the dissipating fields (everything else 
being equal). However, for lower cr, the inflow velocity 
into a given X-point is also reduced, as oc ^fa. Evidently, 
for constant Xcji^ I \f^-, the two opposite effects compen- 
sate, resulting in a fixed low-energy spectral cutoff. 

4.3. Dependence on the Stripe-Averaged Field 

In this section, we explore the dependence of the shock 
properties on the strength of the stripe-averaged field 
{By)\. The magnitude Bq of the field in the cold wind 
is fixed by the value of the magnetization a. A non-zero 
stripe-averaged field is established if the portion of each 
stripe with By = +Bo has a different width than the 

^^ Strictly speaking, the comparison in Fig.JTlis performed at the 
same uipt/y^, namely at the same time in units of the post-shock 
plasma frequency (~ ujp/^/a). However, we stress once again that 
the downstream spectrum does not appreciably evolve in time. 



tates that injection into the SDA process should be easier 
for higher a (at fixed A). 

In summary, we find that both the physics of shock- 
driven reconnection, which ultimately governs the shape 
of the particle spectrum, and the efficiency of SDA are 
controlled by the same parameter Xc/u /V'^j namely the 
stripe wavelength measured in units of the post-shock 
skin depth (~ ^/a c/ujp). We expect that the main prop- 
erties of the shock, and in particular the high-energy end 
of the downstream particle spectrum, should be relatively 
insensitive to variations in A or cr, provided that the ratio 
Xc/uj I \f^ is kept constant. 

In Fig. [7] we investigate this prediction by comparing 
the downstream spectrum among flows that have the 
same value of A^/i^ j \fo — 160/-\/50 ~ 22.5. We vary 
the wind magnetization from cr = 10 up to cr = 100, 
or equivalently the stripe wavelength from A = 72 c/ojp 
up to A = 226c/cijp. The spectra in Fig. It] are shifted 
along the x-axis by 50/cr, to facilitate comparison with 
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region with By = —Bq. If A+ is the length of the former 
and A_ of the latter (with the constraint A+ + A_ = A), 
then we have {By)\/Bo = (A+ — A_)/A. Alternatively, by 
defining a = (A+ — A_)/A+, we have {By)x/BQ — a/(2 — 
\a\). In pulsar winds, \a\ — >■ 1 at latitudes approaching 
the inclination angle between magnetic and rotational 
axes of the pulsar (see the upper panel of Fig. [IJ. The 
special case a = 0, appropriate for the w ind nndplane, 
will be discussed in more detail in 34.3.11 below. 



CT=10 A=320 c/ 



a-0.75 cj„t=3300 



The results presented in the previous sections referred 
to shocks with a moderate value of a = 0.1. As a rep- 
resentative case of high-a flows. Fig. [8] shows the shock 
structure for a = 0.75, corresponding to {By)x/Bo ~ 0.6. 
The stripe wavelength A = 320 c/wp and the wind magne- 
tization cr = 10 are chose n to facilitate comparison with 
the results discussed in §3.1| As in the case a = 0.1, 
the pre-shock flow interacts at first with a fast MHD 
shock (at X ~ 2170 c/wp in Fig. [sl indicated with a ver- 
tical dotted blue line), that triggers magnetic reconnec- 
tion. With increasing a (and so, with increasing {By)\), 
a smaller fraction of the incoming Foynting flux is in 
the form of alternating fields, and as such available for 
dissipation. For a > (which implies A_|_ > A_), recon- 
nection behind the fast shock will eventually annihilate 
the field where By — —Bq, but a finite region where 
the field stays By ~ +Bq will survive. Dissipation of 
the alternating component is completed when the size 
of reconnection islands approaches A_ (or, in general, 
the smallest between A+ and A_), as shown by the 2D 
plots of density and magnetic energy in Fig. ^a) and 
(b), respectively. At this point, a second shock forms (at 
X ~ 2050 c/wp in Fig. Isl) , which would correspond in the 
regime of small a ( i.e., a < 0.1) to the hydrodynamic 
shock discussed in §3.1| For larger a, the speed of this 
second shock increases, from /3sh — 1/3, as appropriate 
for an unmagnetized post-shock plasma, up to /3sh — > 1, 
the shock velocity in a cr ^ 1 flow. In the limit |a| -^ 1, 
the shock that in the regime a < 0.1 was identified as 
"hydrodynamic" becomes eventually degenerate with the 
fast MHD shock. In the case a = 0.75 shown in Fig. [8 
the fast shock stays at a distance of ~ 100 c/uip ahead of 
the main shock, and they both move at /3sh — 0.6. 

Downstream from the main shock {x < 2050 c/wp), 
the residual magnetic flux that has not been annihilated 
on the way from the fast shock survives within stripes of 
high magnetic field, separated by regions of relatively low 
magnetic content (see the 2D and ID plots of magnetic 
energy in Fig. Sth) and (d), respectively; and the ID 
profile of By, black line in Fig. Isle)). The plasma inside 
the high-field stripes is relatively cold, not having gained 
much energy from field dissipation. In contrast, the low- 
field regions are populated with hot particles energized 
by the reconnection field, as shown in the ID profile of 
mean kinetic energy per particle of Fig. ^^) and in the 
longitudinal positron phase space of Fig^j^g). We find 
that the post-shock flow retains a striped structure, with 
magnetically-dominated regions separated by relatively 
wide sheets of hot plasma, for values of |a| > 0.3 (or 
equivalently, for \{By)\\/Bo > 0.15). As \a\ increases, 
the longitudinal extent of the high-field stripes grows, 
at the expense of the regions of hot dense plasma, and 
magnetic pressure dominates over particle pressure in the 
energetics of the downstream fiuid. 

The striped structure of the downstream flow leaves 
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Fig. 8. — Internal structure of the flow at ujpt = 3300, for stripe 
wavelength A = 320c/ajp and magnetization a = 10, zooming 
in on a region around the shock. The stripe-averaged field is 
{By)\/Bo ~ 0.6, corresponding to a = 0.75. The hydrodynamic 
shock is located at a; ~ 2050c/Ldp, and the location of the fast 
MHD shock (x ~ 2170c/Ljp) is indicated with a vertical dotted 
blue line. The following quantities are plotted: (a) 2D plot in the 
simulation plane of the particle number density, in units of the up- 
stream value, with contours showing the magnetic field lines; (b) 
2D plot of the magnetic energy fraction es = B'^ /Sn-fomricoc^ ', 
(c) y-averaged particle number density; (d) y-averaged magnetic 
energy fraction; (e) y-averaged magnetic field By (black line) and 
electric field E^ (red line), normalized to the upstream magnetic 
field Bo; (f) mean kinetic energy per particle, in units of the bulk 
energy at injection; (g) x — 7/3x positron phase space. 



a clear imprint on the post-shock particle spectrum, as 
Fig. |9] shows. There, we present how the particle distri- 
bution changes when the parameter a is varied between 
a = 0.0 and a = 0.95, for fixed wavelength A = 320 c/wp, 
magnetization a = 10, and bulk Lorentz factor 70 = 15. 

For a < 0.3, the average particle Lorentz factor 
downstream from the hydrodynamic shock approaches 
il) llo — cr -I- 1 ~ 11, suggesting full dissipation of mag- 
netic fields (black line in the subpanel of Fig. [9]). The 
shape of the spectrum is independent of ct, in the regime 
a < 0.1. The spectra for a = 0.0 (red curve) and a — 0.1 
(black curve) differ only at the upper end, where the case 
a = 0.1 shows the high-energy bump of shock-drift accel- 
erated particles. This component disappears for a — 0.0, 
when the stripe-averaged field {Ez)x — (3Q{By)x power- 
ing the SDA process vanishes. 

For a > 0.3, the post-shock particle spectrum consists 
of two components. The low-energy peak comes from 
cold plasma residing in the high-field stripes described 
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Fig. 9. — Downstream particle spectrum at LUpt = 3000 for differ- 
ent values of the stripe-averaged field (By) x (or equivalently, of the 
parameter a), in a flow with A = 320c/a;p and a = fO. We vary 
a from o = 0, when the regions with ISy = +Bo and By = — _Bo 
have the same width A+ = A_ = A/2, to the limit a —^ 1, when 
A_ <^ A_|- ~ A. The limit of an unstriped wind is shown for ref- 
erence as a purple line. The black line in the subpanel shows the 
average downstream Lorentz factor as a function of a (with a = 1.0 
referring to the unstriped wind). 



above, whereas the high-energy part is populated by hot 
particles that gained energy from field dissipation. As a 
increases, the fraction of upstream Poynting flux avail- 
able for dissipation decreases, which explains why the 
high-energy component in the spectra of Fig. [9] gets de- 
populated, at the expense of the low-energy part. As 
a consequence, the average particle energy downstream 
from the shock monotonically drops from (7) ~ 7o(cr+l), 
corresponding to the case of full dissipation, down to 
(7) ~ 7o, as expected for negligible field annihilation 
(black line in the subpanel of Fig. [9]). In the limit \a\ -^ 1 
(yellow curve for a = 0.95), most of the downstream 
plasma is contained within the high-field regions, and 
the low-energy component of the spectrum approaches 
the result expected for an unstriped wind (purple line) p^ 
With increasing a, the high-energy component of tne 
downstream spectrum recedes to lower energies (see the 
trend in Fig. |9J from the blue line for a = .45, to the 
yellow line for a = 0.95). As described in [4.1 the maxi- 
mum energy of particles accelerated at X-pomts depends 
primarily on the amount of time available for accelera- 
tion, which is set by the characteristic distance between 
a given X-point and the closest island (in the same cur- 
rent sheet). In turn, this scales with the maximum size of 
reconnection islands, which is determined by the small- 
est between A+ and A_ (here A_ < A+, since a > 0). 
For larger a at fixed A, the extent of A_ decreases, so 
the maximum energy to which particles can be acceler- 
ated by the reconnection electric field decreases. This 
explains the trend in the upper cutoff of the spectra in 



^^ We remark that the low-energy peak in the spectrum of high- 
er flows, including the limit of unstriped wind, is not compatible 
with a Maxwellian (yellow and purple lines in Fig. ^. Here, the 
particle distribution is a ring in momentum space, lying in a plane 
orthogonal to the field. The location of the low-energy spectral 
peak scales linearly with 70, in agreement with the jump conditions 
for a highly magnetized fluid. 



Fig. |9] Given that the location of the high-energy cutoff 
is set by the value of A_ (or equivalently, by the maxi- 
mum size of reconnection islands), one would expect that 
simulations with the same A_, but different choices of A 
and a > 0, should yield similar spectra, at the high- 
energy end. We have tested this, and we find that flows 
with the same A_ show magnetic islands of similar sizes, 
and their downstream spectra almost overlap at high en- 
ergies. Of course, this scaling does not hold for the low- 
energy spectral component, since it is populated by par- 
ticles that were initially in the region with By = +i3oi 
and are therefore insensitive to the value of A_ . 

4.3.1. The Case of Zero Stripe- Averaged Field 

The special case {By)\ — (or equivalently, a = 0), 
which is relevant for the equatorial plane of pulsar winds, 
is worth a deeper investigation. As anticipated in §4.1[ 
hot downstream particles can propagate back from the 
shock into the upstream, and cross at least one wave- 
length of the striped wind, provided that Xc/uj I \f^ ^ 
2 ^2, where ^2 = 1 — 10 is a numerical factor that depends 
weakly on A and g. For wind parameters that satisfy such 
criterion, the motion of hot particles moving into the up- 
stream (the "returning" particles, from now on) will be 
regulated, to zeroth order, by the stripe-averaged fields 
(Bj,)a and (^,)a = /3o(Sj,)a. For a > 0.03 (or equiv- 
alently, {By)\/BQ > 0.015), the returning particles will 
be quickly advected downstream by the stripe-averaged 
fields, after being accelerated via the SDA mechanism 
discussed in fj4.1| On the other hand, for a < 0.03 they 



can propagate as if the upstream flow were unmagne- 
tized. The special case a — 0.0 is taken here as a repre- 
sentative example of the shock structure for a < 0.03. 

Fig. [To] compares the structure of the flow between the 
cases a ~ 0.1 (left panels) and a — 0.0 (right pan- 
els), for the same wavelength A = 80c/wp and mag- 
netization a = 10. The positron phase space of the 
case a = 0.0 (panel (c), right column) shows a dif- 
fuse population of hot particles moving back from the 
shock into the upstream {x > 1050c/cjp), which were 
absent in the case a = 0.1 (panel (c), left column). For 
both a — 0.0 and a = 0.1, as the flow propagates from 
the fast shock (at x ~ 1050c/ajp) to the hydrodynamic 
shock (at X ~ 1000 c/wp), magnetic reconnection trans- 
fers energy from the field to the particles, which in the 
downstream region end up w ith an average Lorentz fac- 
tor (7) ~ 7o(cr + l) (see ^4.1 ). In other words, the physics 



of magnetic reconnection is not sensitive to the value of 
a, in the regime a < 0.1. However, it is only in the 
case a = 0.0 (or a < 0.03, generally speaking), when the 
upstream stripe-averaged fields are small, that the hot 
particles resulting from shock-driven reconnection can 
stream for large distances ahead of the shock, as observed 
in the phase space of Fig. [lO|c) (right panel) . 

When a sufficient number of returning particles builds 
up in the pre-shock region, their counter-streaming with 
respect to the incoming flo w can trigger the filamenta- 



tion (or Weibel) instability (Weibel' 1959| [Medvedev fc 
Loeb|1999| iGruzinov & Waxman 1995|7in tact, the fil- 



amentary structures seen ahead of the shock in the 2D 
plot of Bz (panel (d) , right column) are generated by the 
Weibel instability. A similar pattern appears in the 2D 
plots of electric field Ey (= —/SqBz), particle density, and 
magnetic energy (not shown). In contrast, no signatures 
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Fig. 10. — Comparison of the shock internal structure, for a flow with A = 80c/uip and a = 10, between the case a = 0.1 (left panels) 
and the case a = 0.0 (right panels, corresponding to {By)x = 0). The comparison is performed at uipt = 3000. In both cases, the shock is 
located at a; ~ 1000 c/ujp. As a function of the longitudinal coordinate x, the following quantities are plotted: (a) y-averaged profile of the 
particle density n, in units of the pre-shock value; (b) y-averaged profile of the magnetic energy fraction es; (c) x — 7/3x positron phase 
space; (d) 2D plot on the simulation plane of the out-of-plane component B^ of the magnetic field. 



of filamentation are seen for a = 0.1 (panel (d), left col- 
umn), since in this case no particles can propagate far 
enough ahead of the shock to seed the Weibel instability. 

For both a — 0.0 and a = 0.1, the hydrodynamic shock 
satisfies the jump conditions appropriate for an unmag- 
netized 3D relativistic fluid. However, the fllamentation 
instability appreciably changes the structure of the flow 
for a = 0.0, especially at late times. The density profile 
in the transition region (panel (a)) is much smoother in 
the case a = 0.0 (right column) than it is for a = 0.1 (left 
column). More importantly, the magnetic fields gener- 
ated ahead of the shock by the Weibel instability survive 
in the downstream region of the a = 0.0 shock, as shown 
by the ID profile of magnetic energy in panel (b) (right 
column, at a; < lOOOc/cjp). Instead, for a = 0.1, the only 
magnetic field left behind the shock comes from compres- 
sion of the pre-shock stripe-averaged field, which results 
in a much lower level of post-shock magnetic energy (see 
the left panel of Fig. ^h)). 

The shape of the downstream energy spectrum is also 
significant ly d ifferent between the cases a = 0.0 and a = 
0.1. Fig. [n] follows the time evolution of the particle 
distribution behind a shock with a = 0.0, compared with 
the case a = 0.1. As further discussed in Appendix [Bj 
the downstream particle spectrum in a fiow with a > 0.03 
does not appreciably change with time, so the dotted 
blue line in Fig. 11, which refers to Wpt = 3000 for the 
case a = 0.1, is a good representation of the downstream 
energy spectrum at all times (for a — 0.1). In contrast, 
as shown by the solid lines, the evolution of the particle 
spectrum in the case a — 0.0 is much more dramatic. At 
early times (black line for ujpt — 1500), the spectrum for 
a = 0.0 is remarkably similar to the (time-independent) 
particle distribution of the case a = 0.1 (dotted blue 
line) . This confirms once more that the physics of shock- 
driven reconnection, which is ultimately responsible for 
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Fig. 11. — Time evolution of the downstream spectrum for a 
shock propagating in a flow with A = SOc/iUp, u = 10, and o = 
(i.e., zero stripe-averaged field), from ujpt = 1500 up to uipt = 6000. 
For comparison, the downstream spectrum for a wind with the 
same parameters except for a = 0.1 is shown as a dotted blue line 
for ujpt = 3000 (but it does not appreciably evolve in time). In the 
subpancl, time evolution of the maximum particle Lorentz factor, 
for a = 0.0 (solid line) and a = 0.1 (dotted line). 



shaping the post-shock spectrum at early times, does not 
appreciably depend on a, in the regime a < 0.1. As 
mentioned above, the minor high-energy bump seen in 
the spectrum of a = 0.1 (dotted blue line), but not for 
a — 0.0 (black line), results from particles accelerated at 
the hydrodynamic shock via the SDA mechanism. 

Starting from Wpi = 3000 (solid blue line), the spec- 
trum for a = 0.0 changes significantly. The high-energy 
part of the broad peak that was present at earlier times 
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(at 7 ~ 300) gets de-populated, at the expense of a 
prominent tail of nonthermal particles (at 7 > 10'^) . The 
tail contains approximately a few percent of particles and 
a few percent of the total particle energy. With time, the 
slope of the nonthermal tail becomes flatter (for com- 
parison, the dashed line in Fig. 11 has p = 2.5), 
its high-energy cutoff extends to 



and 



Tigher Lorentz factors, 
as tracked by the solid line in the subpanel of Fig. [TT| 
In contrast, no appreciable evolution in the value of 
the maximum particle Lorentz factor is observed for the 
case a = 0.1 (dotted line in the subpanel). The parti- 
cles in the nonthermal tail seen for a = 0.0 have been 
accelerated via a kind of first-order Fermi mechanism 
(or diffusive shock acceleration, DSA), by scattering off 
the turbulence generated by the Weibel instability (see 
Fig. 10 ^d), right column). Their energization proceeds 



in the same way as for a strictly u nmagnetized shock, 
which was discussed by Spitkovsky ( 2008J . By compar- 
ing the black line (for w^t = 1500) and the blue line 
(for LOpt — 3000) at 7 > 10^, we infer that injection into 
the Fermi process is due primarily to the so-called "ther- 
mal leakage," i.e., the hottest particles from the down- 
stream region can escape ahead of the shock, where they 
take part in the Fermi acceleration cycles. Of course, 
once the fields produced by the Weibel instability grow 
stronger than the background field B^ (see the scale in 
Fig. [lOFd), right column), the Weibel-generated turbu- 
lence m the shock layer may mediate further injection of 
particles into the Fermi process. 

Finally, we remark that, based on the condition 
'^c/wp/v^ ^2^2 discussed above, injection into the 
Fermi process is expected to be more favorable for short 
stripe wavelengths and highly magnetized winds. In 
fact, we observe that the efficiency of Fermi accelera- 
tion, meaning the fraction of particles participating in 
the Fermi cycles, is higher for smaller A (at fixed a) 
and for larger a (at fixed A). This explains why, for 
A = 320c/wp (so, a wavelength longer than the value 



A = 80 c/wp used here) , Fig. M shows that the spectra 
for a = 0.0 and a = 0.1 are alm ost identical, in contrast 
with the conclusions of Fig. Ill] For A — 320c/wp, fewer 
particles can be injected into the Fermi process, and the 
structure of the fiow at uj-pt = 3000 is practically the 
same for all values of a < 0.1. 

4.4. Dependence on the Bulk Lorentz Factor 

Fig. [12] shows the particle spectrum downstream from 
the hydrodynamic shock, for different bulk Lorentz fac- 
tors of the striped wind (from 70 = 3, in red, up to 
7o — 375, in blue). The spectra are shifted along the x- 
axis by I5/70, to facilitate comparison with the reference 
case 7o = 15 (black curve) discussed in the previous sec- 
tions. Once normalized in this way, the spectra overlap 
almost perfectly, suggesting that the physics of shock- 
driven reconnection is not sensitive to the upstream bulk 
Lorentz factor. The value of 70 enters the definition of 
the relativistic plasma skin depth c/ojp (oc ^/to) and of 
the relativistic Larmor radius (r^ — c/ojp/y/a oc ^/to), 
but no residual dependence on 70 remains if length and 
time scales are normalized with respect to c/wp and u]~^ , 
respectively. More specifically, the conditions that con- 



trol the physics of shock-driven reconnection (see j 4.1 
and Petri fc Lyubarsky 2007, ) depend only on the stripe 
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Fig. 12. — Downstream particle spectrum at uj-pt = 3000 for 
different values of the pre-shock bulk Lorentz factor 70, in a flow 
with A = 640c/a;p, a = 10, and a = 0.1. We vary 70 from 70 = 3 
to 70 = 375. Spectra are shifted along the x-axis by I5/70 to 
facilitate comparison with the reference case 70 = 15. The black 
line in the subpanel shows the average downstream Lorentz factor 
as a function of 70. 



wavelength A and the magnetization a, but they do not 
explicitly depend on the wind Lorentz factor. This ex- 
plains why our results are nearly the same (modulo an 
overall shift in energy scale) across a wide range of 70 . 

5. THREE-DIMENSIONAL SIMULATIONS 

The results presented so far are based on 2D experi- 
ments with magnetic field initialized in the plane of the 
simulations. This is the best configuration to study the 
tearing-mode instability, whose wavevector lies in the 
plane of the field, but it artificially inhibits the growth of 
cross-plane instabilities, which may also be important for 
the structure of the shock. In this section, we study the 
3D physics of shock-driven reconnection, with particu- 
lar emphasis on the process of particle energization. The 
simulation setup parallels very closely what we described 
in ^ with the magnetic field oriented initially along the 
y direction. We adopt our fiducial values for cr = 10 and 
a = 0.1, and we investigate two choices for the stripe 
wavelength, A = SOc/wp and A = 160c/wp. For each 
value of A, the extent of the simulation domain along y 
and z is chosen such that our results are not artificially 
affected by the periodicity of our boundary conditions 
(see Appendix ffl . We adopt Ly — L^ — 50 c/uj^ for 



A = 803/_a;p, 
Fig. 



i. 



lOOc/wp for A = 160c/wp. 



T3| shows the structure of the shock at Wpi — 990, 



for a striped wind of wavelength A — IGOc/wp. In agree- 
ment with the results presented in §3.1| for 2D simu- 
lations, the structure of the flow in Fig. [13] shows the 
presence of two shocks, a fast MHD shoclc located at 
X ~ 630 c/wp and a hydrodynamic shock at a: ~ 250 c/wp. 
At the fast shock, the flow is decelerated and compressed, 
as shown by the increase in magnetic energy in Fig. 13 'b). 
The passage of the fast shock through the incoming cur- 
rent sheets initiates magnetic field annihilation. The dis- 
sipation of alternating fields proceeds as the flow propa- 
gates from the fast to the hydrodynamic shock, and little 
magnetic energy remains downstream from the hydrody- 
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Fig. 13. — Internal structure of the flow at ujpt = 990, from the 3D simulation of a striped wind with A = 160c/ajp, cr = 10, and 
a = 0.1. The magnetic field is initialized in the xy plane, as shown by the black arrows in panel (a). The hydrodynamic shock is 
located at x ~ 250c/a;p, and the fast MHD shock at a: ~ 630c/a;p. The following quantities are plotted: (a) 3D plot of the particle 
number density, in units of the upstream value, with color scale stretched to enhance contrast; (b) 3D plot of the magnetic energy fraction 
eg = i3^/87r7omncoc^, with color scale stretched to enhance contrast. 



namic shock (see Fig. 13 |b) at a; < 250c/a;p), in agree- 
ment with our 2D results. As observed in 2D, the par- 
ticle density downstream from the hydrodynamic shock 
(Fig.flsFa)) and the shock velocity are in agreement with 
MHD jump conditions for a relativistic unmagnetized 
fluid. In fact, most of the post-shock en erg y is in kinetic 
form, rather than electromagnetic formFj 

As the flow propagates from the fast to the hydrody- 
namic shock, the structure of the current sheets is af- 
fected by the growth of two competing modes. In the 
xy plane of the magnetic field, the tearing-mode insta- 
bility breaks the current she et into a sequence of high- 
density islands (see Fig. [l3|^a)), separated by X-points 
where the field lines tear and reconnect. In the xz plane 
orthogonal to the field, the folding of the current sheet 
seen in Fig. 13 is governed by the drift-kink instabil- 
ity, driven by the current of fast-drifting plasmas in a 
thin sheet (Daughton 1998 Zenitani &; Hoshino 2005a 
20071. The 2D geometry eniployed in the previous sec- 



tions was chosen to select the tearing mode, and suppress 

^^ As a side note, we point out that the density filamentation 
seen in Fig. |13[ a) ahead of the fast shock is a transient ef fect of the 
passage of the electromagnetic precursor discussed in j ]3.1| The 
precursor wave, and the resulting density filamentation, do not 
appreciably affect the shock structure at late times. 



the drift-kink mode. In the regime of undriven recon- 
nection, the drift-kink instability is generally thought to 
grow faster t han the tearing mode ( Zenitani & Hoshino 



2005b 2007), but for very thin sheets the opposite hier 
archy is observe d, as shown via 3D PIC simulations by 



Liu et ah] (2011 ). In our experiments of shock-driven re- 
connection, the current sheet just behind the fast shock 
{x ~ 580c/wp) shows the pattern of magnetic islands 
typical of the tearing mode, whereas the drift-kink mode 
appears only farther downstream [x < 500c/wp). 

The nature of the instability that dominates the pro- 
cess of field annihilation leaves an imprint on the result- 
ing particle distribution. Zenitani & Hoshino (2007) have 
shown with PIC simulations that held dissipation due to 
the drift-kink instability does not result in nonthermal 
particle acceleration, the plasma just being heated. This 
is because the field lines remain straight, so that all parti- 
cles gain energy at the same rate. Nonthermal particles 



are produced on ly by the tearing mode (e.g., Zenitani 



& Hoshino 20011. In our 2D simulations with in-plane 
fields, i.e., the geometry required to capture the tear- 
ing mode, we have shown that particles are accelerated 
to nonthermal energies by the reconnection electric field. 
They initially move along the z direction parallel to the 
reconnection field, and then they drift along y from a 
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guide field component, though, is no t usually expected in 
the context of pulsar winds (but see Petri fc Kirk|2005). 
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Fig. 14. — In a striped wind with A = SOc/iUp, u = 10, and 
a = 0.1, comparison of the downstream particle spectrum for dif- 
ferent magnetic field configurations: 3D simulation (red) and 2D 
simulations with either in-plane (blue) or out-of-plane (green) mag- 
netic fields. The black line in the subpanel shows the average 
downstream Lorentz factor (7) for the different cases. 



given X-point into the closest magnetic island. In 3D, 
one could argue that the folding of the current sheet in- 
troduced by the drift-kink instability in the xz plane (not 
resolved in 2D) may deflect the accelerating particles out 
of the current sheet, thus suppressing their energization. 

In Fig. |14[ we address this important issue, by com- 
paring for A = 80 c/ujp the post-shock particle spectrum 
of a 3D simulation (red line) to the results of 2D experi- 
ments with in-plane or out-of-plane magnetic fields (blue 
and green lines, respectively). The excellent agreement 
between 3D and 2D in-plane results (red and blue line, 
respectively) suggests that the physics of particle acceler- 
ation by shock-driven reconnection is captured extremely 
well by our 2D simulations with in-plane fieldsrj In turn, 
this implies that, at least for A — 80c/ujp, the maximum 
energy of accelerated particles is constrained by the dis- 
tance between X-points and islands (in the xy plane), 
rather than by the current sheet folding in the xz plane. 
Just upstream of the hydrodynamic shock, the typical 
spacing between X-points and islands scales linearly with 
the stripe wavelength. A similar trend may be expected 
in the characteristic scale of the drift-kink mode, for the 
current sheet just ahead of the hydrodynamic shock (see 
Fig. 13 'b) at a; ~ 350c/wp). If this is the case, our con- 
clusions will also hold for longer stripe wavelengths, as 
we have directly verified in the case A = WQc/utp. A de- 
tailed 3D study of particle acceleration by shock-driven 
reconnection will be presented elsewhere. 

Finally, we point out that the tension force of a guide 
field (along z, in our geometry) can easily stabilize the 
drift-k ink mode, with little or n o effect on the tearing 



mode (Zenitani & Hoshino |2008 ). In the presence of a 
guide field, the 3D physics ot shock-driven reconnection 
should be described very accurately by 2D simulations 
with alternating fields lying in the simulation plane. A 

^* The agreement between red and blue lines in the high-energy 
bump at 7 > 2000 also suggests that the physics of SDA is correctly 
described by our 2D simulations with in-plane fields. 



6. SUMMARY AND DISCUSSION 

We have explored by means of 2D and 3D PIC sim- 
ulations the internal structure and acceleration proper- 
ties of relativistic shocks that propagate in an electron- 
positron striped wind, i.e., a flow consisting of stripes of 
alternating field polarity separated by current sheets of 
hot plasma . This work extends the study of Sironi fc] 
Spitkovskyl (2009), that investigated relativistic shocks 



in a pair plasma carrying uniform magnetic fields. 

We find that a fast MHD shock propagates into the 
pre-shock striped flow, compressing the incoming cur- 
rent sheets and initiating the process of driven magnetic 
reconnection^ via the tearing-mode instability. Recon- 
nection islands seeded by the passage of the fast shock 
grow and coalesce, while magnetic energy is dissipated at 
X-points located in between each pair of islands. When 
reconnection islands grow so big to occupy the entire 
region between neighboring current sheets, the striped 
structure of the flow is erased, and a hydrodynamic shock 
forms. Downstream from the shock, the average particle 
energy is larger than in the pre-shock flow by a factor of 
~ cr, where a is the wind magnetization. In other words, 
the energy stored in the alternating fields has been en- 
tirely transferred to the particles via shock-driven mag- 
netic reconnection, and the downstream fluid behaves 
like an unmagnetized plasma. The only field component 
surviving in the downstream region comes from shock- 
compression of the stripe-averaged pre-shock field (i.e., 
the non-alternating component), if non-zero. In our 2D 
and 3D simulations, we find that complete annihilation 
of the alternating fields occurs irrespective of the stripe 
wavelength A or the magnetiz ation cr, in contrast with th e 
results of the ID analysis by Petri & Lyubarsky (2007). 
In Appendix [X] we show that multi-dimensional studies 
are of paramount importance to correctly capture the 
physics of shock-driven reconnection. 

The main agent of particle energization is the recon- 
nection electric field, as particles drift from a given X- 
point into the closest island. Whether all particles have 
comparable energy gains, or only a few particles are ac- 
celerated to high energies, and the majority stay cold, 
depends sensitively on the properties of the wind. We 
have explored the dependence of the downstream particle 
spectrum on the stripe wavelength A, the wind magneti- 
zation a (or equivalently, the field strength Bo), and the 
stripe-averaged magnetic field {By)\ (or equivalently, the 
parameter a, with a = for {By)\ — 0, and \a\ — >■ 1 in 
the limit \{By)\\ — > Bq). For fixed a = 0.1, we find that 
the shape of the post-shock spectrum depends primarily 
on the combination A/(rx,cr), where r^^ is the relativistic 
Larmor radius in the striped wind. This is just the stripe 
wavelength measured in units of the post-shock plasma 
skin depth. For small values of X/{rL(j) {^ a few tens), 
the spectrum resembles a Maxwellian distribution with 
mean energy (7) ~ 70 cr, where 70 is the bulk Lorentz 
factor of the pre-shock flow. In the limit of very large 
values of A/(ricr) (> a few hundreds), the spectrum ap- 
proaches a broad power-law tail, with a flat spectral slope 
p ~ 1.5. In this regime, the tail extends from 7min ~ 70 



up to 7n 



70 cr 



l/(2-p) 



In terms of the flow structure. 



the former case (i.e., a Maxwellian- like spectrum) is real- 
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ized when all particles pass in the vicinity of an X-point, 
thus gaining energy from field dissipation. In contrast, 
the latter case occurs when most of the particles remain 
far from X-points, thus staying cold, and the energy of 
annihilating fields is transferred to only a few particles, 
which end up dominating the energy content of the flow. 
As |a| increases from 0.1 up to 1, the fraction of up- 
stream Poynting flux in the form of alternating fields be- 
comes smaller. As the amount of field energy available for 
dissipation decreases, the average particle Lorentz fac- 
tor downstream from the shock recedes from (7) ~ "focr, 
as expected for complete field annihilation, down to 
(7) ~ 7o, the result of an unstriped wind. In the op- 
posite limit |a| < 0.1, i.e., in the case of nearly symmet- 
ric stripes, the particles accelerated by the reconnection 
electric field to the highest energies can escape ahead of 
the shock, unaffected by the stripe-averaged field that 
would tend to advect them back downstream. Their 
counter-streaming with respect to the incoming flow can 
seed the filamentation (or Weibel) instability, as it hap- 
pens for shocks propag ating in unmagnetized plasmas 
(Spitkovsky 2005 20081. In turn, the turbulence gener- 
ated by the Weibei instability can mediate particle ac- 
celeration to even higher energies, via a Fermi-like diffu- 
sive process. The nonthermal tail produced by the Fermi 
mechanism is steeper than that given by magnetic recon- 
nection, with characteristic slopes p > 2.5. Typically, a 
few percent of particles are injected into the Fermi pro- 
cess, and the injection efficiency scales in inverse propor- 
tion to the ratio X/{r]^(j). With time, the nonthermal 
tail of Fermi-accelerated particles extends to higher en- 
ergies. For \a\ > 0.1, the propagation of particles far 
ahead of the shock is inhibited by the stripe-averaged 
field, and Fermi acceleration is suppressed. In this case, a 
few percent of particles can get accelerated by the stripe- 
averaged electric field {E;^)\ while gyrating around the 
shoc k front, a process known as shock-drift acceleration 



(i.e., Begelman & Kirk [1990). The tail resulting from 
shock-drift acceleration has a limited energy extent, and 
it does not evolve in time. 

We have confirmed our findings with an extensive set of 
convergence tests. Most importantly, as we have shown 
in 35J 2D simulations with in-plane magnetic fields - the 
setup that we have primarily adopted in this work - can 
reproduce very accurately the shock physics and particle 
spectrum observed in 3D simulations. In other words, 
the current sheet folding caused by the drift-kin k insta- 
bility in the p lane orthogonal to the field (Zenitani & 



Hoshino|2007 ) does not represent a serious constramt tor 



the maximum particle energy achieved in the course of 
the reconnection process. Rather, the maximum energy 
of accelerated particles is limited by the time to drift 
from a given X-point into the closest island. 

These findings can place important constraints on cur- 
rent models of Pulsar Wind Nebulae (PWNe), bubbles of 
synchrotron-emitting plasma powered by the relativistic 
wind of young pulsars. If the magnetic and rotation axes 
of the pulsar are not aligned, the wind consists of stripes 
of opposite magnetic field polarity, alternating with the 
pulsar period. The wind wavelength will be A = 2'!tR]^Ci 
where Rlc = c/il is the so-called light-cylinder radius, 
and fi is the angular frequency of the pulsar. The stripe- 
averaged field vanishes along the equatorial plane (i.e., 
a = there), and the striped structure persists up to a 



latitude equal to the inclination angle between the pul- 
sar's magnetic and rotational axes. At higher latitudes, 
the wind carries a uniform (i.e., non-alternating) field, so 
the physics di scussed in the present work does not apply 
. r^^^, : '' • -^^^^ ' "1992] 



see, instead, Gallant et al. 1992; Hoshino et al. 



Amato fc Arons^ 2006; Sironi fc Spitkovsky ^_2009^ 2011) 



(Jur results suggest that efficient dissipation of the al- 
ternating fields should occur at the termination shock of 
pulsar winds, irrespective of the properties of the flow. 
It follows that current observational constraints on the 



upstream magnetization parameter (iRees & Gunn||1974 



Kennel fc Coroniti||1984a| [Del Zanna et al.||2004[ ) s'Wulc 
be attributed not to the total Poynting flux, but to the 
Poynting flux associated with the stripe-averaged mag- 
netic field, the only component surviving downstream 
from the shock. The upstream flow may be Poynting- 
dominated, provided that most of the Poynting flux is 
in the form of an alternating fleld, which is annihilated 
at the shock. Thus, our results provide a solution to 
the so-called "sigma problem," i.e., that pulsar winds 
should be magnetically-dominated at the light-cylinder 
radius, but kinetically-dominated downstream from the 
termination shock. The idea of shock-driven reconnec- 
tion as a so lution to the sigm a problem was originally 
proposed by Lyubarsky ( 2003 ) . Our work demonstrates 
that complete annihilation of the alternating fields occurs 
regardless of the properties of the wind, in the regime of 
relativistic magnetically-dominated flows. 

In addition, the particle energy spectra extracted from 
our simulations can be used directly to interpret the ra- 
diative signature of PWNe. The radio spectrum of the 
Crab Nebula, the prototype of the class of PWNe, re- 
quires a population of nonthermal particles with a flat 
spectral slope {p ~ 1.5), extending at least across three 
decades in energy, from 10^ MeV up to lO'^ MeV. The 
particle spectrum should be steeper at higher energies, 
with p ~ 2.5, to explain the optical and X-ray flux. 

Our results suggest that, for a slope p ~ 1.5, the width 
of the electron power-law tail resulting from reconnec- 
tion will be ~ a^. To produce a particle distribution 
extending over three decades in energy, the wind mag- 
netization at the termination shock needs to be cr > 30. 
Most importantly, broad particle spectra are produced by 
shock-driven reconnection only if X/(rLa) > a few tens. 
For pulsar winds, the particle number density scales with 
distance from the pulsar as n — nLc{RLc/R)^, where 
the density at the light-cylinder radius is usually writ- 
ten as riLc — kugj- H ere, ngj = O.Blc /2iTec i s the 
Goldreich- Julian density ( Goldreich fc Juliaii||1969 ), and 
K is the so-called multiplicity. Conservation of energy 
along the flow streamlines implies that at the termination 
shock 7o(l-l-cr)K = a;Lc/2i^, where lulc = eBLc/rnc. It 
follows that at the termination shock [R = Rts) 



A 



Attk 



Rlc 
Rts 



(3) 



The constraint X/lr^a) > a few tens required to pro- 
duce broad particle spectra is extremely challenging for 
current models of PWNe and pulsar magnetospheres. 
The radius of the Crab termination shock in the equa- 
torial plane is inferred from X-rav observ ations to be 
Rts ~ 0.1 pc ~ 5 X IO^Rlc ( jHester et alT||2002|). Most 
available models estimate k ~ 10* — 10'^ ( Bucciantini 
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et al.|20TT I , depending on whether or not radio-emitting 
eiectrons are included in the analysis. Based on our find- 
ings, the resulting value of X/irLcr) < 0.01 would yield 
a Maxwellian-like spectrum, at odds with the wide flat 
spectrum required by observations. 

If radio-emitting electrons are produced in the equa- 
torial plane by shock-driven reconnection, a revision of 
the existing estimates of k is required. In this respect, 
we point out that the values of k quoted above are aver- 
ages over latitude, and one cannot exclude that particle 
injection into the pulsar wind is highly anisotropic, with 
multiplicity as large as 10* along the equatorial plane. 
Of course, any given choice of k will constrain 70 and 
a via 7o(l + cf)k = lolc/'^^, a relation that holds at 
each latitude. From the equatorial plane, radio-emitting 
electrons will quickly fill the whole nebula, transported 
by the strong fluid motions observed in MHD models 
of PWNe downstream from the termination shock (e.g., 
Komissarov fc Lyubarsky] [20041 IDel Zanna et al.||2004[ 
Camus et al. 2009 K This would provide a natural ex- 
planation tor the la ck of gradients in the radio sp ectral 
slope of the Crab ([Bietenholz fc Kronberg||1992[ ). Al- 
ternatively, the radio part of the spectrum may be pro- 
duced at higher latitudes, where the termination shock 
is closer to the pulsar, and the ratio in eq. (p| becomes 
larger. However, in this case one should also account for 
the fact that, at high latitudes, the fraction of upstream 
Poynting flux available for dissipation is lower (since \a\ 
is larger), which results in a narrower spectrum, every- 
thing else being fixedp^ 

On the other hand, tKe small value of A/(rLcr) expected 
in the equatorial plane on the basis of current models of 
the Crab Nebula is the most favorable setup for par- 
ticle acceleration via the Fermi diffusive process, with 
efficiency of a few percent by number. In this respect, 
one could interpret the optical and X-ray signature of the 
Crab, which requires a particle spectrum with p ~ 2.5, as 
synchrotron emission from such Fermi-accelerated par- 
ticles. Based on our results, optical and X-ray emit- 
ting electrons should be produced close to the equatorial 



plane of the wind, where |a| < 0.1. If the same equatorial 
wedge is responsible for accelerating the radio-emitting 
electrons, our findings predict that they should outnum- 
ber the optical and X-ray emitting particles by a factor 
of few hundred. This is in agreement with current mod- 
els of PWNe spectra, which req uire k c:^ 10^ ii the radi o 
band is included in the analysis ( Bucciantini et al.|2011 ), 
a value which is two orders of magnitude larger than if 
only higher frequency ba nds are considered (k ~ 10**, 
Kennel fc Coroniti|[l984b| . 

Although our work focuses primarily on the physics 
of PWNe, shocks propagating into striped flows may be 
present also in b lazar jets and gamma-ray bursts (e.g., 
Thompson||2006 ) , where they may provide an appealing 
explanation, based on strong physical grounds, for any 
flat particle distribution inferred from observations. The 
broad-band emission of hotspots in radio galaxies, usu- 
ally interpreted as the termination shocks of relativistic 
jets, indicates that the spectra of accelerated electrons 
need to be fiat (1 < p < 2) below GeV energies (Stawarz 



et al.l2007). In the case of luminous blazar sources, very 



flat electron spectra below Ge V energies are in ferred di- 
rectly from X-ray o bservations ( Celotti fc Ghisellini_2008] 



Sikora et al.||2009[ ). The steeper particle spectrum re- 
quired at ;^GeV energies, with slope p ~ 2.5, could in- 
stead result from Fermi acceleration across the shock, if 
the stripe- averaged field satisfies \a\ < 0.1. 
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APPENDIX 
(A) DEPENDENCE ON THE TRANSVERSE SIZE OF THE SIMULATION BOX 



In Fig. 15 we investigate the dependence of our findings on the transverse size Ly of the computational domain, 
for A = 65lJc/wp, a — 10, and a — 0.1. The agreement between the black curve {Ly = 400c/a;p) and the red line 
{Ly — BOOc/wp) suggests that our results become insensitive to the transverse size of the box (i.e., they converge with 
respect to Ly), for boxes larger than 400c/a;p. Generally speaking, we find that a good criterion for convergence is 
Ly > A/2. When this condition is fulfilled, the growth of magnetic islands behind the fast shock is not artificially 
inhibited by the periodicity of our boundaries in the y direction. Magnetic reconnection then proceeds up to the point 
when islands from two neighboring current sheets will merge, which happens when their size is ~ A/2 (for a < 0.1). 

On the other hand, for Ly < 300 c/wp the growth and coalescence of reconnection islands within a given current 
sheet artificially stops when only one island is left in the sheet, and its size approaches Ly. At this point, a second 
shock forms, loca ted behind the fast shock. For Ly > 400c/cL'p, this would correspond to the hydrodynamic shock 
discussed in ^3.1 where the striped structure of the flow is erased, and field energy is entirely transferred to the 
particles. Instead, in the case Ly < 300c/a;p, complete dissipation of the alternating fields is prohibited by the fact 
that the growth of islands is artificially inhibited by the transverse extent of the box. Here, the post-shock fiow 
retains a striped structure, with regions of hot plasma separated by highly magnetized walls of cold particles. With 
decreasing Ly, less magnetic energy is transferred to the particles behind the fast shock, and the post-shock flow gets 
more dominated by magnetic field energy, with respect to particle kinetic energy (see the de crease in average particle 



energy with decreasing Ly, in the subpanel of Fig. 15 1. This also explains why in Fig. 15 the normalization of the 



low-energy peak, populated by the cold particles in the high-field regions, grows with decreasing Ly, at the expense of 



Also, the termination shock at high latitudes is oblique with 



respect to the wind velocity, a configuration not studied here. 
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Fig. 15. — Downstream particle spectrum at uipt 
with A = 640c/aJn. c r = 10, and o = 0.1. We vary Ly from j 
I&: LyubarskyI l|2007| . The black line in the subpanel shows 



: SOOc/tjp down to Ly = 12.5 c/uj 



in a flow 



^» _ _ _ 

the average downstream Lorentz factor {7) 



p, which approaches the ID setup o f|Petri| 
;or (7) as a function of the box widtti Ly. 



the high-energy component of hot particles that gain ed en ergy from field dissipation. 

As discussed in the main body of the paper (see ^4.1), the maximum size of magnetic islands (which here is set 
by the smallest between A/2 and Ly) correlates with the highest energy at which particles can be accelerated by the 
reconnection electric field. The trend in the upper spectral cutoff of Fig.[T5]can then be explained by the decrease with 
Ly in the maximum size of reconnection islands. In this respect, the spectra in Fig. 1151 interpreted as a monotonic 
trend in the size of reconnection islands, can be directly compared to the spatial sequence of spectra in Fig. ^f ), which 
follow the growth of magnetic islands as the flow propagates from the fast to the hydrodynamic shock. A similar trend 
as in Fig. 15 is also observed in the large-a limit of Fig. l9] (yet for a different wavelength, A = 320c/ajp as opposed to 
A = 640 c/ujp). However, in that case the maximum size of magnetic islands was set by a physical constraint, namely 
the smallest between A+ (where By = +Bq) and A_ (where By = —Bq), whereas here it is the transverse size of the 
box which artificially limits the growth of magnetic islands for Ly < SOOc/wp. 

As the downstream region gets more magnetized with decreasing Ly, the main shock, which would be moving at 
/3sh ^ 1/3 for Lj^ ^ 

The wind parameters employed 
10). In this case, the 



jy > 400c/wp, propagates at a faster veloci ty, eve nt ually catch ing up with the fast MHD shock, in the 
limit Ly ^ A. In this regime, we recover the ID results by Petri & Lyubarsky ( |2007[ ). The 
in this section, A = 640 c/wp and a = 10, satisfy the condition Ac/ui Z''' ;i ^4i/o (where d = 6 - 



ID model by Petri & Lyubarsky (2007) would predict negligible field dissipation, in agreement with our results for a 
very narrow box (yellow hne tor Ly = 12.5 c/ujp] see also the subpanel, for small values of Ly). However, as clarified by 
Fig. |15[ this is just an artificial consequence of the reduced dimensionality of their model, which cannot correctly capture 
the development of the tearing-mode instability, and the resulting growth and coalescence of magnetic islands. In a 
more realistic 2D scenario, complete field dissipation is achieved by the time the flow enters the hydrodynamic shock. 
This stresses that multi-dimensional simulations are essential for our understanding of shock-driven reconnection. 

(B) TIME EVOLUTION OF THE SHOCK 

In Fig. [TBI we follow the time evolution of the shock, from the earliest stages until it reaches a steady state. We 
choose ourfiducial values for A = 640c/wp, a = 10, and (-By) a — 0.05 (corresponding to a = 0.1), but our main 
conclusions hold basically for the whole parameter space we have investigated. At very early times (uppermost row, 
for ujpt = 750), only the fast MHD shock is present (its location at different times is shown as a vertical dotted blue 
line in Fig. 16). The incoming flow decelerates at the fast shock, forming a ring in momentum space lying in the 
xz plane orthogonal to the magnetic field. Since the spread in y-momentum is much smaller than in the other two 
directions, the fluid behind the fast shock can be treated as a 2D relativistic plasma. In fact, the increase in density, 
magnetic energy, and transverse magnetic field across the fast shock, as well as the drop in average kinetic energy per 
particle, are in good agreement with the MHD jump conditions for a 2D relativistic fluid. 

Behind the fast shock, the energy content of the flow at utpt = 750 is dominated by magnetic fields. However, even 
at such early times, we see that current sheets are getting wider (left panel) and populated with hotter particles (see 
the mean kinetic energy per particle, black line in the right panel), as we proceed farther downstream from the fast 
shock. As described in fpj a larger and larger fraction of the cold wind is being channeled into reconnection islands, 
seeded by the passage of the fast shock through the upstream current sheets. At X-points located in between each pair 
of neighboring islands, magnetic energy is being transferred to the particles, which explains why the peaks in mean 
kinetic energy per particle are getting higher with distance behind the fast shock (black line in the right panel), at 
the expense of the mean magnetic energy per particle (red line). Yet, the striped structure, with magnetized regions 
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of cold plasma separated by hot current sheets, still persists downstream from the fast shock. 
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Fig. 16. — Time evolution of the internal structure of a shock propagating in a striped flow with A = GAOc/up, a = 10, and a = 0.1. We 
follow the evolution of the shock from uipt = 750 up to ujpt = 5250 (from top to bottom). Left panels show the jz-averaged particle number 
density, in units of the upstream value. Right panels show the mean kinetic energy per particle (black lines) and the mean magnetic energy 
per particle (red lines). In all panels, the location of the fast MHD shock is indicated as a vertical dotted blue line. 

The shock structure changes significantly at ujpt = 2250 (second row from the top). Behind the fast shock, which 
is now much weaker (see the modest density jump at x ~ 1850 c/wp), reconnection islands have enough time to grow 
and fill the entire region in between neighboring current sheets. At a: ~ TOOc/wp, the striped structure is erased, and 
a hydrodynamic shock forms (see the density profile in the left panel). The density plateau behind the hydrodynamic 
shock is in agreement with the expected jump conditions for a relativistic unmagnetized 3D plasma, namely nd/uu — 4. 
In fact, the energy balance presented in the right panel shows that, by the time the flow crosses the hydrodynamic 
shock, most of the field energy (red line) has been converted to particle kinetic energy (black line). 

A similar picture holds for later times {ujpt — 3750 in the third row from the top, ujpt = 5250 in the lowermost 
row). The density behind the hydrodynamic shock stabilizes at n^ ~ 4nu, and the transfer of energy from magnetic 
to kinetic form becomes more and more complete. The change of adiabatic index due to the transition from a 
magnetically-dominated fluid to a kinetically-dominated plasma results in a lower pressure to drive the fast shock, 
that decelerates (compare the locations of the vertical dotted blue lines at different times) . Between Wpt = 3750 and 
LOpt = 5250, the overall structure of the transition region reaches a steady state. The fast shock propagates ahead of 
the hydrodynamic shock, roughly at a fixed distance (i.e., the speed of the fast shock approaches the velocity f3sh — 1/3 
of the hydrodynamic shock) . The separation between the two shocks obeys the requirement that reconnection islands 
seeded by the fast shock should fill the entire region in between current sheets, by the time the flow arrives at the 
hydrodynamic shock. When the flow structure has reached the steady state, we flnd that the particle spectrum 
downstream from the hydrodynamic shock does not significantly change with time. 
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